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ABSTRACT 

We present synthetic spectral energy distributions (SEDs) for single-age, single- 
metallicity stellar populations (SSPs) covering the full optical spectral range at mod- 
erately high resolution (FWHM=2.3A). These SEDs constitute our base models, as 
they combine scaled-solar isochrones with a empirical stellar spectral library (MILES), 
which follows the chemical evolution pattern of the solar neighbourhood. The models 
rely as much as possible on empirical ingredients, not just on the stellar spectra, but 
also on extensive photometric libraries, which are used to determine the transforma- 
tions from the theoretical parameters of the isochrones to observational quantities. 
The unprecedent stellar parameter coverage of the MILES stellar library allowed us 
to safely extend our optical SSP SEDs predictions from intermediate- to very-old age 
regimes, and the metallicity coverage of the SSPs from super-solar to [M/H] = —2.3. 
SSPs with such low metallicities are particularly useful for globular cluster studies. We 
have computed SSP SEDs for a suite of IMF shapes and slopes. We provide a quan- 
titative analysis of the dependence of the synthesized SSP SEDs on the (in)complete 
coverage of the stellar parameter space in the input library that not only shows that 
our models are of higher quality than other works, but also in which range of SSP 
parameters our models arc reliable. The SSP SEDs are a useful tool to perform the 
stellar populations analysis in a very flexible manner. Observed spectra can be studied 
by means of full spectrum fitting or by using line indices. For the latter we propose 
a new Line Index System (named LIS) to avoid the intrinsic uncertainties associated 
with the popular Lick/IDS system and provide more appropriate, uniform, spectral 
resolution. Apart from constant resolution as function of wavelength the system is also 
based on flux-calibrated spectra. Data can be analyzed at three different resolutions: 
5 A, 8.4 A and 14 A (FWHM), which arc appropriate for studying globular cluster, 
low and intermediate-mass galaxies, and massive galaxies, respectively. Furthermore 
we provide polynomials to transform current Lick/IDS line index measurements to 
the new system. We provide line-index tables in the new system for various popular 
samples of Galactic globular clusters and galaxies. We apply the models to various 
stellar clusters and galaxies with high-quality spectra, for which independent studies 
are available, obtaining excellent results. Finally we designed a web page from which 
not only these models and stellar libraries can be downloaded but which also provides 
a suite of on-line tools to facilitate the handling and transformation of the spectra. 

Key words: galaxies: abundances - galaxies: elliptical and lenticular, cD - galaxies: 
stellar content - globular clusters: general 
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To obtain information about unresolved stellar populations 
in galaxies, one can look at colours, ranging from the UV 



to the infrared. Such a method is popular especially when 
studying faint objects. Since colours are strongly affected 
by dust extinction, one prefers to use spectra when study- 
ing brighter objects. These have the additional advantage 
that abundances of various elements can be studied, from 
the strength of absorption lines (e.g. Rose 1985; Sanchez- 
Bl'azquez et al. 2003; Carretero et al. 2004). Ahhough the 
UV and the near IR are promising wavelength regions to 
study, the large majority of spectral studies is done in the 
optical, mainly because of our better understanding of the 
region. To obtain information, one generally compares model 
predictions with measurements of strong absorption line 
strengths. This method, which is insensitive to the effects of 
dust extinction (e.g., MacArthur 2005) very often uses the 
Lick/IDS system of indices, which comprises definitions of 25 
absorption line strengths in the optical spectral range, origi- 
nally defined on a low resolution stellar library (FWHM>8- 
11 A) (Gorgas et al. 1993, Worthey et al. 1994, Worthey & 
Ottaviani 1997). There is an extensive list of publications in 
which this method has been applied, mainly for early-type 
galaxies (see for example the compilation provided in Trager 
et al. 1998). 

Predictions for the line-strength indices of the integrated 
light of stellar clusters and galaxies are obtained with stellar 
population synthesis models. These models calculate galaxy 
indices by adding the contributions of all possible stars, in 
proportions prescribed by stellar evolution models. The re- 
quired indices of all such stars are obtained from observed 
stellar libraries, but since the spectra in these libraries were 
often noisy, and did not contain all types of stars present in 
galaxies, people used the so-called fitting functions, which 
relate measured line-strength indices to the atmospheric 
parameters {Teg, loggr, [Fe/H]) of library stars. The most 
widely used fitting functions are those computed on the basis 
of the Lick/IDS stellar library (Burstein et al. 1984; Gorgas 
et al. 1993; Worthey et al. 1994) There are, however, alterna- 
tive fitting functions in the optical range (e.g., Buzzoni 1995; 
Gorgas et al. 1999; Tantalo, Chiosi & Piovan 2007; Schiavon 
2007) and in other spectral ranges (e.g., Cenarro et al. 2002; 
Marmol-Queralto et al. 2008). 

The most important parameters that can be obtained 
using stellar population synthesis are age (as a proxy for the 
star formation history) and metallicity (the average metal 
content). There are several methods to perform stellar pop- 
ulation analysis based on these indices. By far the most pop- 
ular method is to build key diagnostic model grids by plot- 
ting an age-sensitive indicator (e.g. H/3) versus a metallicity- 
sensitive indicator (e.g., Mg6, (Fe)), to estimate the age and 
metallicity (e.g. Trager et al. 2000; Kuntschner et al. 2006). 
There are other methods that, e.g., employ as many Lick 
indices as possible and simultaneous fit them in a sense 
(e.g., Vazdekis et al. 1997; Proctor et al. 2004). An alter- 
native method is the use of Principal Component Analysis 
(e.g. Covino et al. 1995; Wild et al. 2009). 

In the last decade the appearance of a generation of 
stellar population models that predict full spectral energy 
distributions (SEDs) at moderately high spectral resolution 
has provided new means for improving the stellar popula- 
tion analysis (e.g., Vazdekis 1999, hereafter V99; Bruzual 
& Chariot 2003; Le Borgne et al. 2004). These models em- 
ploy newly developed extensive empirical stellar spectral li- 
braries with fiux-calibrated spectral response and good at- 



mospheric parameter coverage. Among the most popular 
stellar libraries are those of Jones (1999), CaT (Cenarro et 
al. 2001a), ELODIE (Prugniel & Soubiran 2001), STELIB 
(Le Borgne et al. 2003), INDO-US (Valdes et al. 2004), and 
MILES (Sanchez-Blazquez et al. 2006d, hereafter Paper I). 
Models that employ such libraries are, e.g., Vazdekis (1999), 
Vazdekis et al. (2003) (hereafter V03), Bruzual & Chariot 
(2003), Le Borgne et al. (2004). Alternatively, theoretical 
stellar libraries at high spectral resolutions have also been 
developed for this purpose (e.g.. Murphy & Meiksin 2004; 
Zwitter et al. 2004; Rodn'guez-Merino et al. 2005; Munari et 
al. 2005; Coelho et al. 2005; Martins & Coelho 2007). Mod- 
els that use such libraries are, e.g., Schiavon et al. (2000), 
Gonzalez-Delgado et al. (2005), Coelho et al. (2007). 

There are important limitations inherent to the method 
that prevent us from easily disentangling relevant stel- 
lar population parameters. The most important one is 
the well known age/metallicity degeneracy, which not only 
affects colours but also absorption line-strength indices 
(e.g. Worthey 1994). This degeneracy is partly due to the 
isochrones and partly due to the fact that line-strength in- 
dices change with both age and metallicity. The effects of 
this degeneracy are stronger when low resolution indices are 
used, as the metallicity lines appear blended. Alternative 
indices that were thought to work at higher spectral resolu- 
tions, such as those of Rose (1985, 1994), have been proposed 
to alleviate the problem. 

New indices with greater abilities to lift the 
age/metallicity degeneracy have been proposed with the aid 
of the new fuU-SED models (e.g., Vazdekis & Arimoto 1999; 
Bruzual & Chariot 2003; Cervantes & Vazdekis 2009). These 
models allow us to analyze the whole information contained 
in the observed spectrum at once. In fact, there is a grow- 
ing body of full spectrum fitting algorithms that are being 
developed for constraining and recovering in part the Star 
Formation History (e.g., Panter et al. 2003; Cid Fernandes et 
al. 2005; Ocvirk et al. 2006ab; Koleva et al. 2008). Further- 
more the use of these SSP SEDs, which have sufficiently high 
spectral resolution, has been shown to significantly improve 
the analysis of galaxy kinematics both in the optical (e.g., 
Sarzi et al. 2006) and in the near-IR (e.g.. Falcon- Barroso et 
al. 2003), for both absorption and emission lines. 

Here we present single-age, single-metallicity stellar 
population (SSP) SEDs based on the empirical stellar spec- 
tral library MILES that we presented in Paper I) and Ce- 
narro et al. (2007, hereafter Paper II). These models repre- 
sent an extension of the V99 SEDs to the full optical spec- 
tral range. These MILES-models are meant to provide better 
predictions in the optical range for intermediate- and old- 
stellar populations. In Paper I we provided all relevant de- 
tails for MILES, which has been obtained at the 2.5 m Issac 
Newton Telescope (INT) at the Observatorio del Roque de 
Los Muchachos, La Palma. The library is composed of 985 
stars covering the spectral range AA 3540-7410 A at 2.3 A 
(FWHM). MILES stars were specifically selected for pop- 
ulation synthesis modeling. In Paper II we present a ho- 
mogenized compilation of the stellar parameters (Toff, logg, 
[Fe/H]) for the stars of this library. In fact the parameter cov- 
erage of MILES constitutes a significant improvement over 
previous stellar libraries and allows the models to safely ex- 
tend the predictions to intermediate-aged stellar populations 
and to lower and higher metallicities. 



The models that we present here are based on an em- 
pirical library, i.e. observed stellar spectra, and therefore the 
synthesized SEDs are imprinted with the chemical compo- 
sition of the solar neighbourhood, which is the result of the 
star formation history experienced by our Galaxy. As the 
stellar isochrones (Girardi et al. 2000) - i.e. the other main 
ingredient feeding the population synthesis code - are solar- 
scaled, our models are self-consistent, and scaled-solar for 
solar metallicity. In the low metallicity regime, however, our 
models combine scaled-solar isochrones with stellar spec- 
tra that do not show this abundance ratio pattern (e.g., 
Edvardsson et al. 1993; Schiavon 2007). In a second pa- 
per we will present self-consistent models, both scaled-solar 
and a-enhanced, for a range in metallicities. For this pur- 
pose we have used MILES, together with theoretical model 
atmospheres, which are coupled to the appropriate stellar 
isochrones. 



We note, however, that the use of our (base) models, 
which are not truly scaled-solar for low metallicities, does 
not affect in a significant manner the metallicity and age 
estimates obtained from galaxy spectra in this metallicity 
regime (e.g., Michielsen et al. 2008). For the high metallic- 
ity regime, where massive galaxies reside and show enhanced 
[Mg/Fe] ratios, the observed line-strengths for various pop- 
ular indices (e.g., Mg6, Fe5270) are significantly different 
from the scaled-solar predictions. It has been shown that by 
changing the relative element-abundances, and making them 
a-enhanced, both the stellar models (e.g., Salaris, Groenewe- 
gen & Weiss 2000) and the stellar atmospheres (e.g. Tripicco 
& Bell 1995; Korn, Maraston & Thomas 2005; Cohelo et al. 
2007) are affected. However it has also been shown that base 
models, such as the ones we are presenting here, can be used 
to obtain a good proxy for the [Mg/Fe] abundance ratio, if 
appropriate indices are employed for the analysis (e.g. Ya- 
mada et al. 2006). In fact we confirm the results of, e.g., 
Sanchez-Blazquez et al. (2006b), de la Rosa et al. (2007) 
and Michielsen et al. (2008), who obtain a linear relation 
between the proxy for [Mg/Fe], obtained with scaled-solar 
models, and the abundance ratio estimated with the aid of 
models that specifically take into account non-solar element 
mixtures (e.g., Tantalo et al. 1998; Thomas et al. 2003; Lee 
& Worthey 2005; Graves & Schiavon 2008). 

In Section [2] we describe the main model ingredients, 
the details of the stellar spectral library MILES and the 
steps that have been followed for its implementation in the 
stellar population models. In Section [3] we present the new 
single-age, single-metallicity spectral energy distributions, 
SSP SEDs, and provide a quantitative analysis for assessing 
their quality. In Section |4] we show line-strengths indices as 
measured on the newly synthesized SSP spectra and propose 
a new reference system, as an alternative to the Lick/IDS 
system, to work with these indices. In Section [5] we provide 
a discussion about the colours measured on the new SSP 
SEDs. In Section [6] we apply these models to a set of repre- 
sentative stellar clusters and galaxies to illustrate the capa- 
bility of the new models. In Section [7] we present a webpage 
to facilitate the use and handling of our models. Finally in 
Section [8] we provide a summary. 



2 MODELS 

The SSP SEDs presented here represent an extension to the 
full optical spectral range of the V99 models, as updated in 
V03. We briefly summarize here the main ingredients and 
the relevant aspects of this code. 



2.1 Main ingredients 

We use the solar-scaled theoretical isochrones of Girardi et 
al. (2000). A wide range of ages and metallicities are cov- 
ered, including the latest stages of stellar evolution. A syn- 
thetic prescription is used to include the thermally puls- 
ing AGE regime to the point of complete envelope ejection. 
The isochrones are computed for six metallicities Z =0.0004, 
0.001, 0.004, 0.008, 0.019 and 0.03, where 0.019 represents 
the solar value. In addition we include an updated version of 
the models published in (Girardi et al. 1996) for Z =0.0001. 
The latter calculations are now compatible with Girardi et 
al. (2000). The range of initial stellar masses extends from 
0.15 to 7M0. The input physics of these models has been up- 
dated with respect to Bertelli et al. (1994) with an improved 
version of the equation of state, the opacities of Alexander 
& Ferguson (1994) (which results in a RGB that is slightly 
hotter than in Bertelli et al.) and a milder convective over- 
shoot scheme. A helium fraction was adopted according to 
the relation: Y « 0.23 + 2.25Z. For the TP-AGB phase Gi- 
rardi et al. (2000) adopt a simple synthetic prescription that, 
for example, does not take into account the third dredge-up. 
An improved treatment for this difficult stellar evolutionary 
phase has been recently included by Marigo et al. (2008). 
The effects of such improvements become relevant for the 
near-IR spectral range (see e.g., Marigo et al. 2008; Maras- 
ton 2005). 

The theoretical parameters of the isochrones Tcft, logg, 
[Fe/H] are transformed to the observational plane by means 
of empirical relations between colours and stellar parame- 
ters (temperature, metallicity and gravity), instead of us- 
ing theoretical stellar atmospheres. We mostly use the 
metallicity-dependent empirical relations of Alonso, Arribas 
& Marti'nez-Roger (1996, 1999; respectively, for dwarfs and 
giants). Each of these libraries are composed of ~ 500 stars 
and the obtained temperature scales are based on the IR- 
Flux method, i.e. only marginally dependent on model at- 
mospheres. We use the empirical compilation of Lejeune, 
Cuisinier & Buser (1997, 1998) (and references therein) for 
the coolest dwarfs ( TLff ^ 4000 K and giants (Tcff ^ 3500 K) 
for solar metallicity. For these low temperatures we use a 
semi-empirical approach to other metallicities on the ba- 
sis of these relations and the model atmospheres of Bessell 
et al. (1989, 1991) and the library of Fluks et al. (1994). 
The empirical compilation of Lejeune et al. (1997,1998) was 
also used for stars with temperatures above ~8000 K. We 
apply the metal-dependent bolometric corrections given by 
Alonso, Arribas & Martinez- Roger (1995, 1999; respectively, 
for dwarfs and giants). For the Sun we adopt BCq = —0.12 
and a bolometric magnitude of 4.70. 

We have computed our predictions for several IMFs: the 
two power-law IMFs described in Vazdekis et al. (1996) (i.e 
unimodal and bimodal), which are characterized by its slope 
/X as a free parameter, and the multi-part power-law IMFs of 
Kroupa (2001) (i.e. universal and revised). In our notation 



the Salpeter (1955) IMF corresponds to a unimodal IMF 
with ^ = 1.3. Further details of the IMF definitions are 
given in Appendix A of V03. We set the lower and upper 
mass-cutoff of the IMF to 0.1 and 100 Mq, respectively. 

2.2 MILES stellar library 

To compute the SSP SEDs in the full optical spectral range 
we use the MILES library, for which the main character- 
istics are given in Paper I. The stars of the library were 
selected to optimize the stellar parameter coverage required 
for population synthesis modelling. This is, in fact, one of the 
major advantages of this library. Another important feature 
of the library is that the stellar spectra were carefully flux- 
calibrated. To reach this goal, all the stars were observed 
through a wide (6") slit to avoid selective flux losses due 
to the differential refraction effect, in addition to the higher 
resolution setups used to achieve the blue and red parts of 
the stellar spectra. Telluric absorptions that were present in 
the redder part of the spectra were properly removed. We re- 
fer to Paper I for the technical details on how these aspects 
were tackled. 

To prepare this library for its implementation in the 
models we identified those stars whose spectra might not be 
properly representative of a given set of atmospheric param- 
eters. We used the SIMBAD database for identifying the 
spectroscopic binaries in the MILES sample. A number of 
these stars were removed from the original sample as they 
were found non essential as MILES contains sufficient stars 
with similar atmospheric parameters. We also checked for 
those stars with high signal of variability {AV > 0.10 mag), 
according to the Combined General Catalogue of Variable 
Stars of Kholopov et al. (1998) (the electronically-readable 
version provided at CDS). Some stars with signs of vari- 
ability were removed, when no such high variability was ex- 
pected according to their spectral types and atmospheric 
parameters and, at the same time, we were able to iden- 
tify alternative stars in the library with similar parameters 
and no such sign of variability. Various stars were discarded 
because emission lines were detected in their spectra. For 
various technical reasons we also discarded from the orig- 
inal MILES list an additional number of stars: those with 
very low signal-to-noise in the blue part of the spectrum or 
for which the atmospheric parameters were lacking or had 
problems in the continuum. 

Each stellar spectrum of this selected subsample of stars 
was compared to a synthetic spectrum of similar atmospheric 
parameters computed with the algorithm described in V03 
(Appendix B) (see also 12. 3p employing the MILES stellar li- 
brary. To perform this test we first excluded from the list the 
target star for which we wished to synthesize a similar spec- 
trum. We then compared the synthesized stellar spectrum 
with the observed one. This comparison was found to be 
very useful for identifying stellar spectra with possible prob- 
lems and, eventually, discard them. In case of doubt, par- 
ticularly for those regions of parameter space with poorer 
coverage, we compared the observed spectrum of our star 
with the spectra of those stars, with the closest atmospheric 
parameters, which were selected by the algorithm to synthe- 
size our target star. As an example of this method. Fig. [1] 
shows the results for stars of varying temperature, gravity 
and metallicity. 



According to these careful inspections and tests we 
marked 135 stars of the MILES database. While 60 of these 
stars were removed from the original sample, we kept 75 
stars, which were found to be useful for improving the cov- 
erage of certain regions of parameter space. For the latter, 
however, we decreased the weight with which they contribute 
when we synthesize a stellar spectrum for a given set of at- 
mospheric parameters. In practise this is performed by arti- 
ficially decreasing their signal-to-noise, according to the pre- 
scription adopted in l2.3l For the interested reader we provide 
explanatory notes for these stars in an updated table of the 
MILES sample that can be found in |http:/ /miles. iac.es[ 

2.2.1 Stellar atmospheric parameter coverage 

Figure[2] shows the parameter coverage of the selected sub- 
sample of MILES stars for dwarfs and giants (separated at 
logp = 3.0) compared with the original MILES sample. The 
75 stars with decreased weights are also indicated. The pa- 
rameters of MILES plotted here are the result of an exten- 
sive compilation from the literature, homogenized by taking 
as a reference the stars in common with Soubiran, Katz & 
Cayrel (1998), which have very well determined atmospheric 
parameters. We refer the reader to Paper II for an extensive 
description of the method and the adopted parameters. For 
comparison we show the parameter coverage of some popular 
stellar fibraries (STELIB: LeBorgne et al. 2003; Lick/IDS: 
Gorgas et al. 1993, Worthey et al. 1994; INDO-US: Valdes et 
al 2004; and ELODIE 3.1: Prugniel et al. 2007). For ELODIE 
3.1 we use the parameters determined with the TGMET soft- 
ware (see Prugniel & Soubiran 2004 for details) to maximize 
the number of plotted stars. We adopt the mean value for 
those stars with repeated observations. 

Fig. [5] shows that, at solar metallicity, all types of stars 
are well represented in all libraries. At lower metallicities, 
however, MILES shows a much better coverage of the param- 
eter space for both, dwarfs and giants, than other libraries. 
Particularly relevant is the presence of dwarfs with temper- 
atures above Tofi > 6000 7^ and metallicities [M/H] < -0.5. 
With these stars we are in position to overcome a major lim- 
itation of current population synthesis models based on em- 
pirical stellar libraries providing with predictions for metal- 
poor stellar populations in the age range 0.1 - 5Gyr. Fur- 
thermore, there are sufficient stars to compute models for 
metallicities as low as [M/H] = —2.3, particularly for very 
old stellar populations where the temperature of the turnoff 
is lower. 

Another advantage with respect to the other libraries is 
the coverage of metal-rich dwarf and giant stars that allow 
us to safely compute SSP SEDs for [M/H] = +0.2. Also for 
even higher metallicities MILES shows a good stellar param- 
eter coverage and, therefore, it is possible to compute SEDs 
for [M/H] ~ -1-0.4. These predictions will be presented -and 
we will make them available in our website - in a future pa- 
per, where we feed the models with the Teramo isochrones 
(Pietrinferni et al. 2004), that reach these high metallicity 
values. 

These supersolar metallicity predictions are particularly 
relevant for studying massive galaxies and it is the first 
time they can be safely computed, without extrapolating 
the behaviour of the spectral characteristics of the stars at 
higher metallicities as it is done, e.g., in models based on 
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Figure 1. MILES stars for various representative speetral types are plotted in the different panels in black. In each panel we overplot 
the spectrum corresponding to a star with identical atmospheric parameters that we have computed on the basis of the MILES database 
(except for the star itself) following the algorithm described in Section 12.31 The residuals are given in the bottom of each panel, (see the 
text for details). 



the Lick/IDS empirical fitting functions (e.g. Worthey 1994; 
Thomas et al. 2003). Finally, it is worth noting that these 
new model predictions also extend the range of SSP metal- 
licities and ages of our previous model SEDs based on the 
Jones (1999) (V99) and Cenarro et al. (2001a) (V03) stel- 
lar libraries. A quantitative analysis of the quality of the 
models based on the atmospheric parameter coverage of the 
employed libraries is provided in Section 13.21 



2.3 SSP spectral synthesis 

We use the method described in V99 and V03 for comput- 
ing the SSP SEDs on the basis of the selected subsample of 
MILES stars. In short, we integrate the spectra of the stars 
along the isochrone taking into account their number per 
mass bin according to the adopted IMF. For this purpose. 



each requested stellar spectrum is normalized to the cor- 
responding fiux in the V band, following the prescriptions 
adopted in our code, which are based on the photometric 
libraries described in § 12.11 We refer the interested reader to 
the papers previously mentioned for a full description of the 
method. 

The SSP SED, Sx{t, [M/H]), is calculated as follows: 
/■"It 

SA(i,[M/H]) = / Sx{m,t,[M/-{i])N{m,t)x 

Fv{m,t,[M/B])dm, (1) 

where Sx{m, t, [M/H]) is the empirical spectrum, normalized 
in the V band, corresponding to a star of mass m and metal- 
licity [M/H], which is alive at the age assumed for the stel- 
lar population t. N{m, t) is the number of this type of star. 
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Figure 2. The fundamental parameter coverage of the subsamplc of selected MILES stars is shown in the left panels for dwarfs (top) and 
giants (bottom). Stars with decreased weights are plotted with open circles. The parameters of the original MILES sample are plotted 
in the second column of the panels. We show, for comparison, the parameter coverage of various popular stellar libraries (STELIB, Lick, 
INDO-US and ELODIE). For reference we indicate in each panel the solar metallicity (thin dashed vertical line). 



which depends on the adopted IMF. mi and mt are the stars 
with the smallest and largest stellar masses, respectively, 
which are alive in the SSP. The upper mass limit depends 
on the age of the stellar population. Finally, Fv{m, t, [M/H]) 
is its flux in the V band, which comes from transforming the 
theoretical parameters of the isochrones. It is worth noting 
that the latter is performed on the basis of empirical pho- 
tometric stellar libraries, as described in Section [2. II rather 
than relying on theoretical stellar atmospheres, as it is usu- 
ally done in other stellar population synthesis codes. This 
avoid errors coming from theoretical uncertainties as the in- 
correct treatment of convection, turbulence, non-LTE effects, 
incorrect or imcomplete line list, etc (see Worthey & Lee 
2006). 

We obtain the spectrum for each requested star with 
a given atmospheric parameters interporlating the spectra 
of adjacent stars using the algorithm described in V03. The 
code identifies the MILES stars whose parameters are en- 
closed within a given box around the point in parameter 
space (6*0, log go, [M/H]o). When needed, the box is enlarged 
in the appropriate directions until suitable stars are found 



(for example, less and more metal- rich). This is done by di- 
viding the original box in 8 cubes, all with one corner at that 
point. This reduces the errors in case of gaps and asymme- 
tries in the distribution of stars around the point. The larger 
the density of stars around the requested point, the smaller 
the box is. The sizes of the smallest boxes are determined by 
the typical uncertainties in the determination of the param- 
eters (Cenarro et al. 2001b). In each of the boxes, the stars 
are combined taking into account their parameters and the 
signal-to-noise of their spectra. Finally, the combined spec- 
tra in the different boxes are used to obtain a spectra with 
the required atmospheric parameters by weighting for the 
appropiate quantity. For a full description of the algorithm 
we refer the reader to V03 (Appendix B). 

It is worth noting that, within this scheme, a stellar 
spectrum is computed according to the requested atmo- 
spheric parameters, irrespective of the evolutionary stage. 
We are aware that this approach is not fully appropriate for 
the TP-AGB phase, where O-rich, C-rich and stars in the 
superwind phase are present. In fact, in Marmol-Queralto 
et al. (2008) we have found non negligible differences in the 
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Figure 3. Comparison of MILES temperatures with those of Gonzalez-Hernandez & Bonifacio (2009) as a function of temperature (left 
panel) and metallicity (right panel). Giants are represented by red solid circles (81 stars) whereas dwarfs are represented by open blue 
circles (134 stars). The error weighted linear fit to all the stars is plotted in the two panels with a solid line. The dashed line in each panel 
shows the mean error weighted offset obtained from all the stars AT^ff = 51(±10)K, whereas the dotted lines represent the standard 
deviation (140 K) from this offset. 



CO bandhead at 2.3 /im between 19 AGB stars, all from the 
MILES library, and RGB stars of similar parameters. This, 
in principle, could introduce an error in our predictions as 
very few C-rich, O-rich or superwind phase TP-AGB stars 
can be found in MILES. However, the contribution of such 
stars to the total fiux budget in the V-band is just a few 
percent, mainly for stellar populations of intermediate-ages 
(0.1 - 1.5 Gyr; see e.g., Bruzual 2007). 

2.3.1 Effects of systematic variations in the adopted 
stellar parameters 

Percival & Salaris (2009) have recently shown that system- 
atic uncertainties associated with the three fundamental stel- 
lar atmospheric parameters might have a non negligible im- 
pact on the resulting SSP line-strength indices. The inter- 
ested reader is refered to that paper for details and for a 
suite of tests showing the effects of varying these parame- 
ters. 

In particular, a relatively small offset in the effective 
temperature of 50-100 K, which is of the order of the sys- 
tematic errors in the conversion from temperature to colours 
used here (Alonso et al. 1996) may change the age of a 14 Gyr 
stellar population by 2-3 Gyr and alliviate the so-called zero- 
point problem for which the ages of the globular clusters 
(GCs) are older than the most recent estimations of the age 
of the Universe. Furthermore, they show that, in many cases, 
there is a mistmatch in scales between the underlying models 
(the isochrones) and the adopted stellar library in the stellar 
population models. We showed already in Paper II that the 
Alonso et al. (1996, 1999) photometric library did not show 
any offset with the MILES atmospheric parameters. There- 
fore, there is not a missmatch in scales between the models 
and the stellar library in the models presented here. 

However, inspired by this work we decided to compare 
our stellar parameters with those from several recently pub- 
lished works, e.g., Gonzalez-Hernandez & Bonifacio (2009), 
Casagrande et al. (2006), Ramirez & Melendez (2005), to 



see if our parameters should be corrected. We show here, 
for ilustrative purposes, the comparison of our parameters 
with those using the new temperature scale of Gonzalez- 
Hernandez & Bonifacio (2009), whom derived Teff using the 
infrared flux method (Blackwell et al. 1990), as it is the work 
with a larger number of stars in common with MILES (215) 
and because it is the one showing larger differences. 

In Fig.|3] we plot the difference in temperature for 
the dwarf and giant stars in common between Gonzalez- 
Hernandez & Bonifacio (2009) and MILES as a function 
of temperature and metallicity. We performed linear fits 
weighted with the errors obtaining the following relations: 

ATeff = -116(±80) + 0.0312(±0.0148)refiMiLES (2) 
and metallicity 

ATeff = 7(±17) - 36.56(±11.79)[Fe/H] (3) 

obtaining an rms of 138 K in the two cases. The mean tem- 
perature offset of 59 K and 54 K for dwarfs and giants, re- 
spectively. This is in good agreement with the offsets ob- 
tained by Gonzalez-Hernandez & Bonifacio (2009) when 
comparing their temperature scale with that of Alonso et 
al. (1996) and Alonso et al. (1999). Finally we did not find 
any significant offset in metallicity or gravity among these 
samples. 

We computed an alternative SSP model SED library by 
transforming the MILES temperatures to match this scale. 
These models can be used to better assess the uncertainties 
involved in the method. We comment on these models in 
Sections [5l and [6". 1.11 However, we do not find any strong rea- 
son for adopting the Gonzalez-Hernandez & Bonifacio (2009) 
since the accuracy at high metallicities is worse than in the 
Alonso et al. (1996,1999) (the temperature scale of Gonzalez- 
Hernandez & Bonifacio (2009) has been optimized for low 
metallicities while the errors in the temperature determina- 
tion for metal rich giant stars are very large). Furthermore, 



as we already mentioned above, the temperature scale of the 
isochrones and stellar libraries in our models do not show 
any offset. 



3 MILES SSP SEDS 

Table [l] summarizes the spectral properties of the newly syn- 
thesized SSP SEDs. The nominal resolution of the models is 
FWHM=2.3 A, which is almost constant along the spectral 
range. For applications requiring to know the resolution with 
much greater accuracy, we refer the reader to Fig. 4 of Pa- 
per I. This resolution is poorer than that of the models we 
presented in V99 (FWHM=1.8 A), but this is compensated 
by the large spectral range AA 3540.5-7409.6 A, with flux- 
calibrated response and no telluric residuals. For computing 
the SEDs we have adopted a total initial mass of 1 Mq. Ta- 
ble [1] also summarizes the SSP parameters for which our 
predictions can be safely used. We discuss this issue in Sec- 
tion 



3.1 Behaviour of the SSP SEDs 

In Fig. |4]we show the new spectral library of SSP models for 
different ages and metallicities. All the spectra are plotted at 
the nominal resolution of the models (FWHM=2.3A). The 
SSP age-sequence shows weaker 4000 A break and stronger 
Balmer line strengths as the age decreases. The largest 
Balmer line strengths, however, are reached for ages of ~0.4 
Gyr. On the other hand the metallic absorption lines be- 
come more prominent with increasing age. It can be seen 
that the continuum of the models with larger ages is heav- 
ily lowered by the strengthening of the metallic features 
(line blanketing). This is the reason why the Lick/IDS in- 
dices were defined using psedocontinua instead of real con- 
tinua (Worthey et al. 1994). In the SSP SED metallicity- 
sequence it can be seen that Balmer line-strengths decrease, 
and metallicity-sensitive features get stronger for increasing 
metallicity. These two SSP SED sequences clearly illustrate 
the effects of the age/metallicity degeneracy. 

In Fig. [S] we illustrate the effects of varying the IMF 
slope (fi) and the IMF shape. The left panel shows a se- 
quence of SSP spectra of solar metallicity and 10 Gyr, where 
the slope is varied for the Unimodal IMF. For old stellar 
populations this IMF type shows the largest variations as 
the fraction of low mass stars vary dramatically when vary- 
ing the IMF slope. We see that the reddest part of the spec- 
trum shows the largest sensitivity to this parameter. These 
variations do no evolve linearly as a function of /i. In fact 
the SSP SEDs corresponding to the flatter IMFs look rather 
similar. On the other hand, the effects caused by varying 
fi quickly become stronger for the steepest IMFs, notably 
redward of ~5500 A, including the Na line at ~5800A and 
the TiO molecular band around ~6200 A and ~6800A. The 
right panel shows that the effects of varying the IMF slope 
are much less significant for the Bimodal IMF. This is be- 
cause for this IMF type the contribution from stars with 
masses lower than 0.6 M0 is decreased. Since this IMF more 
closely resembles the Kroupa IMFs, it is probably more real- 
istic. Finally, we note that the SSP SEDs for the two Kroupa 
(2001) IMFs and the Bimodal IMF with slope 1.3 (as well 
as the Salpeter IMF) do not show significant differences. 



3.2 Reliability of the SSP SEDs 

In this section we try to provide some quantitative way of 
estimating the reliability and the quality of the synthesized 
SSP SEDs, and to estimate the size of the errors due to an 
incomplete coverage of the stellar parameter space of the 
input library. We perform this evaluation as a function of 
SSP age, metallicity and IMF. To do this, we have taken 
advantage of our algorithm for synthesizing a representative 
stellar spectrum for a given set of atmospheric parameters 
(see Section [2.3^ to compute for each SSP a parameter, Q, 
as follows 
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(4) 

We integrate the term enclosed within the square-brackets 
along the isochrone, from the smallest stellar mass, mi, to the 
most massive star, mt, which is alive at the age of the stellar 
population t, weigthing by the flux in the V band, Fvi, and 
the number of stars, Ni, according to the adopted IMF. A'^^^ 
is the number of stars that have been found within cube j in 
parameter space (where j = 1, . . . , 8, see Section [2. 3[) . For a 
star with given 6, log (7 and [M/H], the size of these cubes is 
given by (pel Aog gi and (/"[m/h];- ^e^, <^iogg^ and o-[m/h]„ is 
the minimum uncertainty in the determination of 9, log g and 
[M/H], respectively. Following V03 we adopt ag^ = 0.009, 
""log 9m ~ ^-^^ ^^'^ '^l^/'nim ~ O-O^- ^hc slzc of the smallest 
cube is obtained by multiplying the minimum uncertainty by 
the factor Xs- Therefore, according to the square-bracketed 
term, for each star along the isochrone, the larger the total 
number of stars found in these 8 cubes and the smaller the 
total volume, the larger the quality of the synthesized stellar 
spectrum is. 

We also compute, for each SSP, the minimum accept- 
able value for this parameter, Qm- Our algorithm ensures 
that at least we find stars in 3 (out of 8) cubes to synthe- 
size a stellar spectrum (see Appendix B of V03 for details). 
Furthermore the max;imum enlargement that is allowed for 
each cube is assumed to be 1/10 of the total range covered 
by the stellar parameters in our library. Therefore the sizes 
for the largest cubes are aoM = 0.17, o-\ogg^ — 0.51 and 
<^IM/H]j^ = 0.41. Note that for erg we also assume the con- 
dition that 60 ^ Tcff ^ 3350 isT, which corresponds to the 
adopted ae limiting values when transforming to the Teg 
scale. Therefore to compute Qm we will assume for each 
stellar spectrum along the isochrone the following three con- 
ditions that must be satisfied: we find stars in three cubes, 
ii) there is a single star for each of these three cubes, which 
implies that Am = Ej=i-^/~3, and m) it is necessary to 
open the size of these cubes, generically, /p x o^y^, fp =0.34, 
i.e. the star is found to be located within 1 cr of o-p^^, con- 
sidering the two directions at either side of the requested 
parametric point. By substituting these values in EqUwe 
obtain 
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Figure 4. Left panel: Solar metallicity SSP spectra of different ages (in Gyr), increasing from top to bottom. Right panel; SSP spectra of 
10 Gyr and different metallicities, increasing from top to bottom. All the SSP spectra are plotted at the nominal resolution of the models 
(i.e. FWHM=2.3 A). All these models are computed adopting a Kroupa Universal IMF. For clarity, the SSP spectra have been shifted by 
arbitrary amounts. 



Qr. 



8/p ^i=l 



NjFy- 



"[M/Hlm , 



(5) 



Then a quality measure of the SSP SED, due to the atmo- 
spheric parameter coverage of the stellar library feeding the 
models, is given by the normalized parameter 



Qn = 7^ 



(6) 



Although the absolute value of Qn depends on the assump- 
tions adopted for computing Qm, the Qn parameter allows 
us to compare the quality of a given SSP SED compared to 
the SEDs computed for other SSP parameters. In addition, 
Qn allows us to compare the quality of SSP SEDs synthe- 
sized on the basis of different stellar spectral libraries, given 
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Figure 5. Left panel: solar metallicity SSP spectra of lOGyr, computed using a Unimodal IMF of increasing slope from top to bottom. 
The Salpeter case, i.e. fi =1.3, has been plotted in the right panel of Fig. (4] The SSP spectra are plotted at the nominal resolution of the 
models (i.e. FWHM=2.3A). Right panel: solar metallicity SSP spectra of lOGyr for different IMF shapes. From top to bottom: Kroupa 
Universal, Kroupa Revised and Bimodal with increasing slopes, as indicated in the panel. All the spectra have been shifted by arbitrary 
amounts for clarity. 



Table 1. Spectral properties and parameter coverage of the synthesized MILES SSP SEDs 



Spectral properties 



Spectral range 

Spectral resolution 

Linear dispersion 

Continuum shape 

Telluric absorption residuals 

Units 



AA 3540.5-7409.6 A 
FWHM = 2.3A, cr = 54km s"! 
O.gA/pix (51.5km s"!) 
Flux-scaled 
Fully cleaned 



Fa/L 



0^ 



^M 







: 3.826 X lO^^erg.s-l 



SSPs parameter coverage 



IMF type 

IMF slope (for unimodal and bimodal) 

Stellar mass range 

Metallicity 

Age ([M/H] = -2.32) 

Age ([M/H] = -1.71) 

Age (-1.31 [M/H] ^ -fO.22) 



Unimodal, Bimodal, Kroupa universal, Kroupa revised 
0.3 - 3.3 
0.1 - 100 Mq 
-2.32, -1.71, -1.31, -0.71, -0.41, 0.0, -1-0.22 
10.0 < t < ISGyr (only for IMF slopes ^ 1.8) 
0.07 < i < 18 Gyr 
0.06 < i < 18 Gyr 



/J. — 0.5 /i,= 1.5 yCi = 2.3 /2 = 3.3 
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Figure 6. The quality parameter, Qn, as a function of the SSP age (in Gyr) for different metaUicities, indicated by the line types as 
quoted in the inset. For computing the SEDs we use a Unimodal IMF with the slope indicated on the top of each panel. SSPs with Q„ 
values larger than 1 can be safely used (see the text for details). 



the fact that the stellar spectra of these libraries are of suffi- 
ciently high quality. Furthermore these prescriptions are the 
same that we use for computing a stellar spectrum of given 
atmospheric parameters, which have been extensively tested 
in V03. 

Figure [§] shows the value of Q„ as a function of the 
SSP age (in Gyr) for different metaUicities (different line 
types) and adopting a Unimodal IMF, with its slope increas- 
ing from the left to the right panel. SSP SEDs with Q„ values 
above 1 can be considered of sufficient quality, and therefore 
they can be safely used. The panels show that, as expected, 
the higher quality is achieved for solar metallicity. For the 
Salpeter IMF, i.e. the second panel, Q„ reaches a value of 
~5 for stellar populations in the range 1-10 Gyr. Outside 
this age range the quality starts to decrease as low-mass and 
hotter MS stars are less numerous in the MILES database. 

A similar behaviour is found in Fig. |S] for the SSPs 
of metaUicities [M/H] = -1-0.22, -0.41 and -0.71, aU with 
Qn values in the range 3-4. Despite the fact that for 
all these metaUicities the value of Q„ sharply decreases to- 
ward younger stellar populations, the computed SSP SEDs 
can be considered safe. For lower metaUicities we obtain 
significantly lower Q„ values, which gradually decrease to- 
ward younger ages. For old stellar populations we obtain 
Q„ values around 2 and around 1.5 for [M/H] — —1.31 
and [M/H] — —1.71, respectively. Interestingly, according to 



our quantitative analysis, the SSP SEDs synthesized for 
[M/H] — —2.32 and ages above 10 Gyr are barely acceptable 
for IMFs with slopes ^ = 0.3 and = 1.3. These SEDs 
are particularly useful for globular cluster studies. Note the 
sharp increase of Q„ observed for ages around 1 Gyr. This 
should be mostly attributed to the fact that for stars of tem- 
peratures larger than 9000 K we do not take into account 
the metallicity of the available stars for computing a stellar 
spectrum of given atmospheric parameters (see V03 for de- 
tails). Therefore, the quality of the SEDs synthesized for the 
younger stellar populations is lower than stated by Q„. 

Fig. |S] also shows that Qn decreases slightly with in- 
creasing IMF slope. In fact, the largest values are obtained 
for the flatter IMF {fj, =0.3), as a result of the smaller contri- 
bution of low mass stars, which are less abundant in MILES. 
This explains why Qn drops significantly for the steepest 
IMF. It is worth noting that the SSP SEDs that are synthe- 
sized for the lowest metallicity cannot be considered safe for 
IMFs steeper than Salpeter. 

For comparison we show in Fig. [7] the Qn values ob- 
tained for two of the most popular libraries that are used in 
stellar population synthesis, i.e. Lick and STELIB. The first 
library is mainly used to predict Lick/IDS line indices via 
the empirical fitting functions of Worthey et al. (1994) (e.g., 
Worthey 1994; Vazdekis et al. 1996; Thomas et al. 2003). 
STELIB is the library feeding the models of Bruzual & Char- 
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Figure 7. Same test as in Fig. |6] but for SSP SEDs computed on 
the basis of the atmospheric parameters corresponding to the Lick 
(left panel) and STELIB (right panel) stellar spectral libraries (see 
the text for details). To perform this test we selected the MILES 
spectra for the stars in common with MILES, whereas for the rest 
of the stars we synthesized their spectra with our interpolating 
algorithm, using the MILES spectral database (see the text for 
details). 

lot (2003). For this purpose we used MILES to assemble a 
stellar spectral library resembling Lick, and another one re- 
sembling STELIB. For each of these libraries we selected 
the MILES spectra for the stars in common with MILES 
with listed stellar parameters (i.e. 213 for Lick and 117 for 
STELIB). For the remaining stars, which are not available in 
MILES, we synthesized their spectra with our interpolating 
algorithm of Section [13] using MILES. 

The two panels of Fig. [7|show that the obtained Qn val- 
ues are significantly lower than the ones obtained for MILES 
for all metallicities and ages. The left panel shows that the 
SSP predictions based on the Lick library are not reliable for 
metallicities lower than ~ —1.3. The right panel shows that 
this limit occurs at higher metallicity (~ —0.7) for any pre- 
diction based on STELIB. OveraU the STELIB SSP models 
are of much lower quality than the predictions based on the 
Lick library, judging from the parameter coverage of these 
libraries. Note that there is a large dispersion in the Q„ val- 
ues obtained for old stellar populations. This feature is less 
prominent in Fig. (6] due to the fact that the distribution of 
stellar parameters is far more homogeneous for MILES. 

Table [T] summarizes the ranges in age, metallicity and 
IMF slopes for which our MILES SEDs can be considered 
safe. We do, however, strongly recommend to inspect Fig. [6] 
for a more accurate description of the quality of the SSP 
SEDs. Although not shown, we find that the Qn values ob- 



tained for the two Kroupa IMF shapes and the Bimodal IMF 
with slope 1.3 are virtually identical to those shown in the 
second panel of Fig. [6] For the Bimodal IMF the Qn val- 
ues also decrease with increasing IMF slope, but the effect 
is much lower than is the case for the Unimodal IMF. 



3.3 Model SED comparisons 

For an extensive comparison of our newly synthesized SSP 
SEDs with the models of other authors we refer the reader 
to the recently published paper of Koleva et al. (2008) and 
Cid-Fernandes & Gonzalez 

(2010fl 

Koleva et al. compared 
three SSP spectral libraries: GALAXEV (Bruzual & Chariot 
2003), Pegase-HR (Le Borgne et al. 2004) and ours. Whereas 
the two former codes use the Padova 1994 stellar evolu- 
tion models (Bertelli et al. 1994), ours employ the Padova 
2000 models (Girardi et al. 2000), whose red giant branch 
are about 50-200 K hotter. Each of these codes use differ- 
ent prescriptions for transforming the theoretical stellar pa- 
rameters to the observational plane. The main difference is 
that our models mostly rely on empirical relations, rather 
than on model atmospheres. However, the key difference be- 
tween these models is the stellar libraries that they use, i.e., 
STELIB, ELODIE (version 3.1) and MILES for GALAXEV, 
PEGASE-HR and ours, respectively. Koleva et al. (2008) do 
this comparison by means of the full spectrum fitting ap- 
proach, using two different algorithms: NBURSTS (Chilin- 
garian et al. 2007), which performs a parametric non-linear 
fit, and STECKMAP (Ocvirk et al. 2006ab), which employs 
a non-parametric formalism. 

Koleva et al. (2008) find that the quality of the mod- 
els is very dependent on the atmospheric parameter cover- 
age of the stellar spectral libraries. The two fitting methods 
show consistent results when our SEDs are compared to the 
PEGASE-HR SEDs, obtaining residuals of the order of 1 per 
cent. As expected, the largest discrepancies are found for low 
metallicities ([M/H] < —0.7), as shown in their Fig. 2. How- 
ever in Fig. [2] we show that MILES provides a significantly 
better coverage than ELODIE 3.1 for this low-metallicity 
regime both, for dwarfs and giants. Furthermore the last two 
panels of Fig. 2 of Koleva et al. show that the discrepancy 
between these two model sets are larger for SSPs with ages 
smaller than ~5 Gyr, which can be explained by the lack of 
hot metal-poor dwarfs in the ELODIE library. The fits per- 
formed by these authors reveals that the worst results are 
obtained when either PEGASE-HR or our models are com- 
pared with the SEDs of Bruzual & Chariot (2003). Koleva et 
al. conclude that this is due to the poor stellar atmospheric 
coverage of STELIB, particularly for the non solar metal- 
licity regime, as is also shown in our Fig. [2] This has been 
recently confirmed by Cid-Fernandes & Gonzalez (2009) by 
comparing the model spectra to Milky Way globular clusters. 

When comparing our MILES models with the models 
that we presented in V99 we show in Fig. [8] a representative 
SSP spectrum of solar metallicity and 10 Gyr, with the cor- 
responding V99 model over plotted. The latter is an updated 
version of the original V99 models, which has been computed 
with the same code used to synthesize the MILES SSP SEDs. 



^ We provided these authors with a preliminary version of our 
MILES models to perform this comparison 
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Figure 8. SSP SEDs of solar mctallicity and 10 Gyr computed 
with the MILES (thin sohd line) and Jones (1999) stellar libraries, 
i.e. the V99 model (thick dotted line). We adopted a Unimodal 
IMF of slope 1.3. The V99 model, which has a higher resolution 
(FWHM=1.8A) than MILES, was smoothed to match the res- 
olution of the MILES model (FWHM=2.3A). The residuals arc 
plotted in the lower panel with the same scale. 



The V99 SSP spectrum, with a resolution of 1.8 A (FWHM), 
was smoothed to match the resolution of the MILES model, 
i.e., 2.3 A (FWHM). The residuals are plotted on a similar 
scale. The residuals show a systematic trend in the two nar- 
row spectral ranges covered by the V99 models. This must 
be attributed to errors in the flux-calibration quality of the 
latter, as was reported in V99. The same pattern is seen in 
the residuals of Fig. [Q] which shows a similar comparison at 
metallicity [M/H] = —1.7. Contrary to Fig. [S] the residuals 
of Fig. [9] show much larger variations at shorter wavelength 
scales. This translates into larger index strength variations 
in this low metallicity regime, showing that the Jones (1999) 
library feeding the V99 models lacks important stars. This 
comparison illustrates how relevant it is to synthesize SSP 
SEDs with stellar libraries with good atmospheric parame- 
ters coverage, such as MILES. 



4 DEFINING A NEW STANDARD SYSTEM 
FOR INDEX MEASUREMENTS 

A major application of the SSP SEDs is to produce key line- 
strength indices that can be compared to the observed val- 
ues. In practise this has been the most popular approach 
for studying in detail the stellar content of galaxies. So far 
the most widely used indices are those of the Lick/IDS sys- 
tem (e.g., Gonzalez 1993; Trager et al. 1998, 2000; Vazdekis 
et al. 1997; j0rgensen 1999; Sanchez-Blazquez et al. 2003, 
2006a,b,c, 2009). The standard Lick/IDS system of indices 
was defined on the basis of a stellar spectral library (Gor- 
gas et al. 1993; Worthey et al. 1994), which contains about 
430 stars in the spectral range AA4000-6200 A. The Lick/IDS 
system has been very popular because the stellar library con- 
tains stars with a fair range of stellar parameters and, there- 
fore, is appropriate to build stellar population models. The 
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Figure 9. Similar comparison as in Fig. |8] but for metallicity 
[M/H] = -1.7. 



most popular application of this library has been the use of 
empirical fitting functions, which relate the index strengths, 
measured at the Lick/IDS resolution, to the stellar atmo- 
sphere parameters. Because the value of the Lick/IDS indices 
depend on the broadening of the lines, authors with spectral 
data who wish to use the population models based on the 
Lick/IDS fitting functions (e.g., Worthey 1994; Vazdekis et 
al. 1996; Thomas et al. 2003), or compare them with previ- 
ously published data, need to transform their measurements 
into the Lick/IDS system. As the resolution of the Lick/IDS 
library (FWHM~8-llA) is much lower than what is avail- 
able with modern spectrographs, the science spectra are usu- 
ally broadened to the lower Lick/IDS resolution. This is per- 
formed by convolving with a Gaussian function, whose width 
varies as a function of wavelength, as the resolution of the 
Lick/IDS system, apart from being low, also suffers from an 
ill-defined wavelength dependence (see Worthey & Ottaviani 
1997). This effect is particularly significant for systems with 
small intrinsic velocity dispersions, such as globular clusters 
and dwarf galaxies, for which part of the valuable informa- 
tion contained in the higher resolution galaxy spectra is lost. 

After that, it is still necessary to correct the measure- 
ments for the line broadening due to the velocity dispersion 
of the stars in the integrated spectra of galaxies. The cor- 
rection should be calculated by finding a good stellar tem- 
plate of the galaxy by combining an appropriate set of stel- 
lar or synthetic spectra which matches the observed spec- 
trum. This is because this correction is highly sensitive to 
the strength of the index. Despite of this, many studies still 
use single polynomial obtained with a single star as a tem- 
plate, that may have very different index-values than those 
on the observed spectra (see Kelson et al. 2006 for a discus- 
sion of the systematic effects associated to this procedure). 

Due to the fact that the shape of the continuum in the 
Lick/IDS stars is not properly calibrated, the use of the Lick 
system requires a conversion of the observational data to 
the instrumental response curve of Lick/IDS dataset (see 
the analysis by Worthey & Ottaviani 1997) that will affect, 
predominantly, the broad indices. This is usually done by ob- 



serving a number of Lick stars with the same instrumental 
configuration as for the science objects. Then, by comparing 
with the tabulated Lick index measurements, the authors 
find empirical correction factors for each index. This experi- 
ment shows that even for the narrow indices, differences exist 
between the Lick/IDS stars and other libraries (see Worthey 
& Ottavianni 1997; Paper I). The origin of these differences 
is not clear, but they need to be corrected. 

Furthermore, the spectra of this library have a low ef- 
fective signal-to-noise ratio due to the significant flat-field 
noise (Dalle Ore et al. 1991; Worthey et al. 1994; Trager et 
al. 1998). This translates into larger random errors in the 
indices, much larger than in present-day galaxy data. It is 
worth noting that the accuracy of the measurements based 
on the Lick system is often limited by these transformations, 
rather than by the quality of the galaxy data. 

4.1 A new Line Index System: LIS 

With the advent of new stellar libraries, such as MILES, 
with better resolution, better signal-to-noise, larger number 
of stars and a very complete coverage of the atmospheric 
parameter space, it is about time to revisit the standard 
spectrophotometric system at which the indices are mea- 
sured. It is not our objective to tell the authors at which 
resolution the indices need to be measured. In fact, we rec- 
ommend to use the flexibility provided by the model SSP 
SEDs and compare the measurements performed on each ob- 
ject with models with the same total broadening. However, 
this might not be a valid option in many cases, e.g., when 
a group of measurements is to be represented together in 
the same index-index diagram. It is also convenient to agree 
on a standard resolution at which to measure line-strengths, 
so that measurements can be compared easily with previous 
studies, with colleagues, and to avoid systematic errors that 
might affect the conclusions. One of the reasons why the 
Lick/IDS system is still popular nowadays is because most 
previous studies have calculated their indices in this sys- 
tem. Although it is a good standard procedure to compare 
new measurements with previous works, the fact that those 
are transformed into the Lick/IDS systems perpetuates the 
problem. For this reason, we propose in this section a new 
Line Index System, hereafter LIS, with three new spectral 
resolutions at which to measure the Lick indices. Note that 
this new system should not be restricted to the Lick set of in- 
dices in a flux calibrated system. In fact, LIS can be used for 
any index in the literature (e.g., for the Rose (1994) indices), 
including newly defined indices (e.g., Cervantes & Vazdekis 
2009). 

We provide conversions to transform the data from the 
Lick/IDS system to LIS as well as tables with index measure- 
ments in the new system for popular samples of Milky Way 
globular clusters, nearby elliptical galaxies and bulges. For 
this purpose we use index tables already published by other 
authors, and transform them into the new system. To cal- 
culate the conversion, we have broadened the MILES stars 
in common with the Lick/IDS library (218 stars in total) to 
the new standard resolutions defined for both, globular clus- 
ters and galaxies, and compare the indices measured in both 
datasets. Third-order polynomials were fitted in all cases. 
These transformations are given in Appendix [A] In a forth- 
coming paper, we will present empirical fitting functions for 
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Figure 10. SSP SED of solar metallicity, 10 Gyr and Kroupa 
Universal IMF, plotted for different resolutions. From top to bot- 
tom, FWHM=2.3, 5.0, 8.4 and 14.0 A. The latter three resolutions 
correspond to the new system of indices proposed here (LIS). 



the Lick indices at these new standard resolutions proposed 
here. A preliminary version of these fitting functions can be 
found in Martfn-Hernandez et al. (2007), which were com- 
puted at the resolution 2.3 A (FWHM). We concentrate here 
on the Lick/IDS suite of indices because it has been the most 
widely used in the literature. 

The definition of this new system using the MILES stel- 
lar library has two main practical advantages: 

• The resolution of the MILES stars is constant as a func- 
tion of wavelength. As a result, the user will not have to 
degrade the spectra with a wavelength dependent function. 

• Additional offsets will not have to be applied to the 
data, as long as the measurements are flux calibrated (in a 
relative sense). Therefore the system is universal, and the 
index measurements straightforward to perform. 

It is very difficult to choose an appropriate resolution at 
which to define a new system, as it depends on the intrin- 
sic broadening of the objects and the instrumental resolu- 
tion. The Lick/IDS indices have been measured in a variety 
of objects, from globular clusters to giant elliptical galax- 
ies, covering a very wide range in velocity dispersion, and 
in data with a variety of instrumental resolutions. The ap- 
propriate resolution also depends on which indices are going 
to be measured; while some Lick/IDS indices do not change 
their sensitivity to age and metallicity with broadening, oth- 
ers (e.g. Ca4227, Fe4531) do. For these reasons, we consider 
it appropriate to define three different standard resolutions. 
The first one has been specifically chosen for stellar clus- 
ters and dwarf galaxy data, while the other two are more 
appropriate for working with massive galaxy spectra. 

4.1.1 LIS-S.OA. 

When choosing a standard resolution at which to compare 
indices of GCs, we were guided by the resolution of exist- 
ing high-quality GC data. Schiavon et al. (2005) presented 
a spectral library of 40 Galactic GCs with a resolution of 



~ 3.1 A(FWHM). This represents an important comparison 
database for extragalactic GC systems, and a logical upper 
limit for our choice of resolution. However, the majority of 
extragalactic spectral data for GCs have lower resolution. 
Data from the Gemini/GMOS collaboration (e.g., Norris et 
al. 2008; Pierce et al. 2006ab) exhibits a typical FWHM 
of ~5A. Very similar to that are the VLT/FORS spectra 
presented by Puzia et al. (2004). Keck/LRIS data from the 
SAGE^I group has been generally taken at slightly higher 
resolutions ~2.8-4A (e.g. Strader et al. 2005; Beasley et al. 
2006; Cenarro et al. 2007; Chomiuk et al. 2008; Beasley et 
al. 2009). Since it is considerably more straightforward to 
degrade higher resolution data to lower resolution than it 
is the attempt to correct lower resolution indices to higher 
resolution, we suggest using a resolution (FWHM) of 5 A in 
GC index comparisons, a choice which encompasses all of the 
above data. This choice reflects a factor of 2 improvement in 
resolution over the Lick/IDS system. 

Authors willing to compare their data with the models 
presented here will only have to broad their spectra to a res- 
olution of 5 A (FWHM) (this corresponds to cr = 127kms~^ 
at 5000 A) and measure the indices in there. If they wish 
to compare with previous published data in the Lick/IDS 
system, they can use the conversion provided in Table \K\ 



4.1.2 LIS-8.4A. 

For galaxies we have decided to define two different stan- 
dard resolutions. The first one, discussed in this subsection, 
is fixed at resolution FWHM=8.4A (this corresponds to 
a = 214kms"^ at 5000A). The reasons for this choice are 
twofold: 

• It matches, roughly, the resolution of the Lick/IDS 
system in the H/3 - Mg triplet region. The indices in this 
regions (e.g., H/3, Mgb, Fe5270) are the most widely used in 
the literature and, therefore, the measurements in the new 
system will not differ considerably from previous measure- 
ments. 

• 200kms~^ is roughly the mean velocity dispersion 
for the early-type galaxies in the Sloan Digital Sky Survey 
(Bernardi et al. 2003a,b,c). 

This system is particularly appropriate for studies 
of dwarf and intermediate-mass galaxies, where the total 
broadening of the spectra (cyf^^^ + Cgai) is not higher than 
a = 214 km s"^ at 5000 A. In this case, the author will only 
have to broad their spectra to the proposed total a using a 
Gaussian broadening function and measure the indices. 

If, on the contrary, the total broadening of the object 
is higher than this value, the author has the option of using 
the system LIS-14.0A, or to correct his/her measurements to 
a total broadening of 8.4 A(FWHM) using polynomials (see, 
e.g., Kuntschner 2000; Sanchez-Blazquez et al. 2006a). 

In a similar way as for LIS-5A, previous data in the 
Lick system can be converted into this new system using the 
polynomials given in Table \K2\ 



4.1.3 LIS-I4.0A. 

The third system is defined at a resolution of FWHM=14 A 
(this corresponds to a =357kms"^ at 5000 A). LIS-14.0A 
has been designed especially for the study of massive galax- 
ies. In the spectra obtained for these galaxies the broadening 
is dominated by the velocity dispersion of the stars. The total 
broadening is higher than the standard Lick/IDS resolution 
for which, in order to compare with the models, usually one 
had to perform a correction due to velocity dispersion of the 
galaxy. As mentioned above, this correction is one of the 
largest sources of systematic errors in the measurement of 
the indices (e.g. Kelson et al. 2006). 

For studies that only analyse the stronger Lick in- 
dices in a sample of galaxies with velocity dispersions 
a > 200 kms~^, we suggest to broaden all the spectra 
by an amount such that the total resolution (crtotai = 
^'^i'nBt + '^gai + <o.i) matches FWHM=14A, before mea- 
suring the indices, where (Tinat represents the instrumental 
resolution and cTgai the velocity dispersion of the galaxy. This 
will help us to avoid to have to perform any further correc- 
tion to the indices due to broadening. 

Fig [To] shows the synthetic spectra of an SSP with 
10 Gyr and solar metallicity broadened to the resolution 
of the different LIS- systems. The polynomials required to 
transform the indices on the Lick/IDS system to this reso- 
lution are provided in Table IXSl 

4.2 Brief notes about system transformation 

Most spectrographs keep the resolution constant in units of 
wavelength, rather than velocity. Therefore, the broadening 
of the spectra, if done in linear wavelength scale, should be 
always performed using gaussians of a fixed number of pixels. 
If the broadening of the data has to be made with a Gaus- 
sian with a constant width in kms~^. then the broadening 
in the spectra should be performed on a spectrum with a 
logarithmic wavelength scale. 

It is common practise in studies of stellar populations 
to observe stars in common with the Lick/IDS library. This 
is done to derive offsets to transform the data onto the spec- 
trophotometric system of the Lick librarjjf]. We stress here 
that this is necessary, only because the Lick/IDS spectra suf- 
fer from a number of calibration problems. However, these 
ofi'sets will not correct properly any uncertainty in the cal- 
ibration of the science spectra, as the position of the line- 
index on the detector changes for extragalactic objects due 
to their redshift. 



^ http:/ /www. ucolick.org/~brodic/sagcs/SAGES_WeIcome. html 



^ For most indices, the corrections are better described by linear 
relations instead of offsets, see appendix in Sanchez-Blazquez et 
al. (2009) 



Table 2. Lick indices measured in the Line Index System for globular clusters and galaxies from selected references. 
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The complete version of this table can be found in the electronic version of the journal. 
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Figure 11. H/3 vs. [MgFe] diagnostic diagram in the LIS-S.OA 
system. Asterisks represent the globular cluster sample of Schi- 
avon et al. (2005), whereas the solid circles represent the clusters 
of Puzia et al. (2002). We used the indices listed in Table [2] (the 
electronic version). For the models we used a Kroupa Universal 
IMF. 



4.3 Reference tables 

One of the reasons many authors keep using the Lick/IDS 
system is to compare with previous studies, which are usu- 
ally presented in this system. We provide here with a set of 
indices, obtained from several references already transformed 
to the new systems. The data are available at the CDS and 
in the electronic version of the paper. 

We have selected several samples of early-type galaxies 
and globular clusters from the literature and transformed 
the line-strength indices to the new systems defined here. 
This is intended to provide a benchmark of indices properly 
transformed into LIS and measured on high quality data that 
other authors can use to compare with. The selected samples 
are the following: 

• Puzia et al. (2002). This sample comprises 12 Galactic 
globular clusters. Lick indices were measured from long-slit 
spectra obtained with the Boiler & Chivens Spectrograph 
mounted on the 1.5m telescope in La Silla. The wavelength 
coverage of the spectra ranges from 3400 to 7300 A with a 
final resolution of ~6.7 A. All the spectra were observed with 
a slit with of 3'.'0. We transformed the Lick indices into the 
new system using the transformations provided in this paper 
(Table mj. 

• Schiavon et al. (2005). The sample consists of 40 galac- 
tic globular clusters, obtained with the Blanco 4-m telescope 
at the Cerro Tololo observatory. The spectra cover the range 
'-3350-6430A with 3.1 A (FHWM) resolution. The signal-to- 
noise ratio of the flux calibrated spectra ranges from 50 to 
240 A"^ at 4000 A and from 125 to 500 A"^ at 5000 A. The 
sample has been carefully selected to contain globular clus- 
ters with a range of ages and metallicities. The indices pro- 
vided in Table [2] have been measured directly in the spectra 
degraded to the LIS-5.0A resolution. 

• Trager et al. (2000). This sample consist of 40 galax- 
ies drawn from that of Gonzalez (1993) intended to cover, 
relatively uniformly, the full range of velocity dispersion. 



line-strength and color displayed by local elliptical galax- 
ies. The line-strength indices (i.e. H/3, Mgb, Fe5270, Fe5335) 
were measured from long-slit data obtained with the CCD 
Cassegrain Spectrograph on the 3m Shane Telescope of Lick 
Observatory. The spectral features were measured within 
re/8 aperture. The indices were transformed into the new 
system applying the transformations from Tables IA2I and 

ESI 

• Sanchez-Blazquez et al. (2006a). This sample consists 
of 98 early-type galaxies including ellipticals and lenticulars, 
spanning a wide range in central velocity dispersion (a). The 
sample contains galaxies in different environments: isolated 
galaxies, galaxies in groups, and in several clusters. The 
central spectra were extracted within an aperture 0.6 kpc 
(assuming Ho=70 kms~^). The spectra were corrected from 
the presence of nebular emission using GANDALF (Sarzi et 
al. 2006), instead of using the [OIII]A-based correcion pub- 
lished in the original paper. We use the models provided 
in this paper as templates. The spectral range covers from 
-3500/3700 (depending on the run) to ~5200/5700 A, with 
spectral resolution ranging from 3.6 to 6.56 A (FWHM) (see 
original paper for more details). The indices provided in Ta- 
ble [2] have been measured in the emission- free spectra previ- 
ously degraded at the different LIS resolutions. For those 
galaxies with a a larger than the LIS resolution an ad- 
ditive correction was made using the kinematic templates 
(the combination of SSP-models that better reproduced the 
galaxy spectra obtained in the measurement of a using ppxf 
(Cappellari & Emsellem 2004)), following the prescriptions 
in Kelson et al. (2006). The correction was obtained as the 
difference in the indices measured in these templates de- 
graded at the LIS-resolution and at the a^^i + af^^ , where 

(Tgai is the velocity dispersion of the galaxy and (Tinst the in- 
strumental resolution. 

• Yamada et al. (2006). The sample consists of 14 galax- 
ies, morphologically classified as E in the Virgo cluster. The 
galaxies where selected with three criteria: (1) they had to 
be classified as elliptical galaxy in the RC3 catalogue (de 
Vaucouleurs et al. 1991); (2) they had to be free from strong 
emission lines; (3) they have to cover a wide range of lumi- 
nosity evenly distributed along the color-magnitude of the 
Virgo cluster. These extremely high quality spectra, which 
have been extracted within an aperture of re/10, have signal- 
to-noise ratio well above >100A~^ for all the galaxies, and 
about -5001"^ for some of them. The spectra cover a wave- 
length range from 3800 to 5800 A, and the resolution ranges 
from 2.0 A to 3.1 A, depending on the observational runs, 
which have been performed at the William Herschel Tele- 
scope (WHT; 4.2m) and Subaru (8.2m) telescopes. For full 
details of the sample we refer the reader to the original pa- 
per. Three galaxies have a velocity dispersion larger than 
the resolution of the LIS-8.4A and the indices in this system 
were corrected in the same way as in the SB06b sample 

• Kuntschner et al. (2006). This sample consists of 48 
E/SO galaxies in the SAURON survey of early-type galaxies 
(de Zeeuw et al. 2002) at the William Herschel 4.2m tele- 
scope in La Palma. We use here the line strength measure- 
ments (i.e. H/3, Fe5015, Mg6) averaged within re/8 aperture. 
The spectral coverage of the observed SAURON integral- 
field observations ranges from 4800 to 5380 A at resolution 
~4.2 A (FWHM). The indices of Table [2] were obtained ap- 
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Figure 12. H/3 vs. [MgFc] and H/3 vs. [MgFcSO] diagnostics diagrams in LIS-8.4A. In tlic left panel wc plot the galaxy samples of Trager 
et al. (2000) (solid circles), Yamada et al. (2006) (stars) and Sanchez-Blazquez et al. (2006a) (open triangles). In the right panel we 
plot the galaxy sample of Kuntschner et al. (2006) (asterisks) and 36 galaxies from the sample of Sanchez-Blazquez et al. (2006a) (open 
triangles), as for all these galaxies the two iron lines required to calculate the [MgFe] index fall outside the covered spectral range (i.e. 
Fe5270 and Fe5335). 
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Figure 13. Same as in Fig.[T2]but for LIS-14.0A. 



plying the transformations of Table [X2l and [X3l t o the indices 
transformed into the Lick systems provided in the original 
paper. 

Table [2] shows a portion of the tables to illustrate its for- 
mat and content. The table contains the Lick/IDS indices 
measured at the different LIS resolutions with the index def- 
initions provided in Trager et al. (1998). The full tables are 
available in the electronic version of the paper. 

Fig. [11] shows the commonly used index-index diagnostic 
diagram H/3 vs. [MgFe] in the new LIS-5.0A system for the 
Puzia et al. (2002) and Schiavon et al. (2005) globular cluster 
samples. We use the indices in Table (2] Note that we do 



not plot the model indices for SSPs younger than 10 Gyr 
for [M/H] — —2.3, according to our Q„ parameter analysis. 
In Fig. [Il and Fig. [13] we show the H/3 vs. [MgFe] and H/3 
vs. [MgFe50] diagrams for the four galaxy samples both in 
LIS-8.4A and LIS-14.0A. Note that for all the galaxies of 
Kuntschner et al. (2006) and for 36 galaxies in the Sanchez- 
Blazquez et al. (2006a) sample it was not possible to plot the 
[MgFe] index, as either one or the two iron indices (Fe5270, 
Fe5335) required to calculate this index were outside the 
covered spectral range. Note that we used GANDALF (Sarzi et 
al. 2006) to remove nebular emission from the galaxy spectra 
of Sanchez-Blazquez et al. (2006a) and Yamada et al. (2006). 
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Figure 15. Time evolution of various Balmer line-strength in- 
dices measured on the SSP SEDs of solar metallicity and Kroupa 
Universal IMF smoothed at different resolutions as in Fig. 1141 
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Figure 17. Variation of Balmer line-strength indiees with metal- 
licity predicted by our single stellar population models of 10 Gyr 
and Kroupa Universal IMF. Different line-styles show predictions 
for indices measured at different spectral resolution, as in Fig. 1141 



We refer the reader to all these original works for a detailed 
discussion of the obtained results. 



4.4 MILES SSP SED index measurements 

Figure O shows the predicted time evolution of our solar 
metallicity and Salpeter IMF models for a selection of line- 
strength indices measured at four different spectral resolu- 
tions: the nominal resolution of the models, i.e. FWHM — 
2.3 A, and the resolutions of LIS-5.0A, LIS-8.4A, and LIS- 
14.0A (i.e., FWHM=5, 8.4 and 14 A, respectively). It can 
be seen that the sensitivity of some indices to variations 
in resolution is barely noticeable. This is the case for the 
wide molecular index definitions, i.e., CN, Mg and TiO. For 
some lines, however, broadening the spectra means loosing 
part of the information contained in them (e.g. Ca4227). The 
most critical situation when analyzing data by means of di- 
agnostic diagrams based on these indices happens when the 
relative variation of the index with age depends on resolu- 
tion. Authors should keep this in mind in order to choose 
the right system and diagnostics to work with the data. Our 
recommendation is to use the most straightforward approach 
possible, which is to work with the model SEDs smoothed 
to match the total resolution of the data. This approach al- 
lows the user to take full advantage of the whole information 
contained in the data (see V99 for an extensive discussion). 
Interestingly we note that the time evolution for most Lick 
indices shows a minimum for intermediate-age stellar pop- 



ulations in the range 0.1-1.0 Gyr, as already mentioned in 
Sec. 13.11 This is particularly evident for the CN indices. 

Fig. US] is similar to Fig. [14] but here several Balmer 
index definitions are shown. We plot the standard Lick H/3 
index and the four higher-order Balmer indices of Worthey & 
Ottaviani (1997). We also include two of the age indicators of 
Vazdekis & Arimoto (1999), the recently defined H/3o index 
of Cervantes & Vazdekis (2009) and the D4000 break (e.g., 
Gorgas et al. 1999). We see that for all the Balmer index 
definitions the dependence on resolution is marginal or vir- 
tually nil (e.g. the 4000 A break). The two U-y^ of Vazdekis 
& Arimoto (1999) are specifically defined to be insensitive to 
resolution. For this reason we only plot, for H7125, the values 
corresponding to the resolution of LIS-5.0A and for H7200 
the ones for LIS-8.4A as these two resolutions fall within the 
stability range of these indices, respectively. We did not in- 
clude the third index defined in Vazdekis & Arimoto (1999), 
H7275, as the resolution of LIS-14.0A is larger than is al- 
lowed for this index. Interestingly, being almost negligible, 
we see that, whereas the sensitivity to resolution for the new 
H/3o index is slightly larger than that for the standard H/3 
Lick index for old stellar populations, the situation is reverse 
for ages smaller than ~lGyr. The higher-order Balmer line 
index definitions of Worthey & Ottaviani (1997) also show 
larger resolution sensitivity for younger stellar populations. 
Note that all Balmer indices, except the D4000 break, peak 
at ~0.3 Gyr and, therefore, the fits obtained on the basis of 
these indices can have two solutions. Kauffmann et al. (2003) 
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Figure 14. Time evolution of several Lick line-strength indices measured on our new SSP SEDs for solar metallicity and Kroupa Universal 
IMF, smoothed to different resolutions. Solid line: resolution of LIS-5.0 A (FWHM = 5 A); dashed line: resolution of LIS-8.4A (FHWM 
= 8.4A); dashcd-dotted lino: resolution of LIS-14.oA (FWHM = 14 A). 



take advantage of the behaviour of these indices to define di- 
agnostics for constraining the Star Formation History. 



Figure [16] and Fig. \T7\ shows the variation of the line- 
strength indices as a function of metallicity predicted by our 
new models for stellar populations of 10 Gyr and Salpeter 
IMF. As in the previous figures the predictions are shown 
for three difTerent resolutions; FWHM= 5.0, 8.4 and 14 A. 
We see a resolution-dependent differential behaviour as func- 
tion of metallicity for several indices, such as for example 
Ca4227, Fe5335, Fe5782. This shows that a lot of care has to 
be taken to match the resolution of data and models when 
comparing both as, otherwise, artificial trends can be de- 
rived (see Kelson et al. 2006). Fig. 1171 shows the somewhat 
opposite dependence on the resolution of two H/3 index defi- 
nitions: H/3o shows larger resolution dependence than H/3 for 
higher metallicities. Note however the insensitivity to metal- 
licity achieved by this H/3o definition for [M/H] > —1, which 
is not seen for the commonly used H/3 Lick index. Such in- 
sensitivity is also shown by the two H7 indices of Vazdekis 
& Arimoto (1999) plotted here. We are not describing here 
in detail the index bahaviour as function of metallicity and 
age, but the metallicity insensitivity of the two TiO molec- 
ular indices for [M/H] > —0.5 is particularly noteworthy. 



4.5 H/3: comparison with the fitting function 
predictions 

In Fig. [18] we focus specifically on the behaviour of the most 
popular age indicator, the Lick H/3 index, as a function of 
age (plotted in linear scale, for ages >1 Gyr) and metallicity. 
In this figure we compare the index computed on the basis 
of the empirical fitting functions of Worthey et al. (1994), 
using the Vazdekis et al. (1996) code, as updated in this 
work and transformed to the LIS-8.4A system following Ta- 
ble |A2l with the H/3 index measured on the SSP SEDs, once 
smoothed to the LIS-8.4A resolution. We see that the models 
based on MILES provide larger values for solar and super- 
solar metallicities and for [M/H] < — 1 than those based on 
the Lick/IDS fitting functions. This implies that when using 
MILES models the obtained ages are older by about ~2 Gyr 
for SSPs around ~10Gyr. 

Another interesting effect at metallicities [M/H] < — 1 
is that, for ages above ~12Gyr, H/3 starts increasing instead 
of decreasing as expected. This effect, not seen at higher 
metallicities, is attributed to the contributions of the Hori- 
zontal Branch (HB) stars (e.g. Maraston & Thomas 2000). 
Therefore the H/J index, rather than indicating the age of the 
stellar population due to their sensitivity to the T^g of the 
turnoff stars, become sensitive to the HB morphology. This 
result is in agreement with Maraston & Thomas (2000). We 
see, however, that this sudden increase in the strength of the 
H/? index is less visible for the new predictions for the three 
lowest metallicities plotted in this figure. In fact, the effect 




virtually disappears for [M/H] < —2. This result would have 
important implications for the interpretation of integrated 
spectra of composite stellar systems, for which the Balmer 
lines, such as the H/3 index, are the primary age diagnostics. 
In fact, Maraston & Thomas (2000) have suggested that a 
an spread in metallicity of an old stellar populations can ac- 
count for most of the dispersion of H/3 values observed in the 
ellipticals family - due to the higher H/3 index of the metal 
poor stars- without the need of a young population as had 
been previously suggested in other studies (e.g. Gonzalez 
1993). As we obtain with MILES higher H/3 strengths, the 
required fraction for the metal-poor component should be 
smaller than is suggested by Maraston & Thomas (2000). 

It is very important to note that in Fig. [18] we are com- 
paring the same synthesis code, with identical HB prescrip- 
tions. The different behaviour at low metallicities between 
the two sets of predictions should be fully attributed, there- 
fore, to the limitations of the Lick stellar library, as it be- 
comes evident when comparing the Qn values obtained for 
these two sets of predictions in this low metallicity regime 
(see the second panel of Fig. [6] and first panel of Fig. [TjiJ. 



^ In fact, guided by the obtained Qn parameter, we only plotted 
the H/3 values for ages larger than 10 Gyr for the MILES models 
of [M/H] = —2.3. However we plotted the H/3 values obtained for 
the Lick/IDS fitting functions for [M/H] = —1.7, despite the fact 
that their corresponding Qn values are lower than 1 (see the first 
panel of Fig. [3 
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Figure 18. Comparison of the H/3 index computed on the basis 
of the empirical fitting functions of Worthey et al. (1994), using 
the models of Vazdekis et al. (1996), as updated in this work 
and transformed to the LIS-8.4A system following Table lA2l with 
the index measured on the SSP SEDs, once smoothed to the LIS- 
8.4A resolution. For the lowest metallicity we only plot the results 
for SSPs older than 10 Gyr according to our qualitative analysis 
shown in the second panel of Fig. (6] However, for the sake of the 
discussion we also show the predictions based on the Lick/IDS 
fitting functions for [M/H] = — 1.7 despite we consider they do 
not have enough quality due to the low number of stars at those 
metallicities (see the text for details). 



5 MILES SSP COLOURS 



One of the main advantages of MILES is that the stellar 
spectra were carefully flux-calibrated. In fact we obtained 
for each star of the library a wide-slit spectrum to avoid 
selective flux losses due to the differential refraction effect 
(see Paper I for details) . This allows us to measure accurate 
colours from the synthesized SSP SEDs. This is performed by 
convolving the SED with the corresponding filter responses. 
To obtain the Johnson B — V broad band colour we use the 
filters of Buser & Kurucz (1978), with the zero points set 
with Vega. To obtain the zero point for the B and V filters, 
we used two different Vega SEDs from Hayes (1985) and 
from Bohlin & Gilliland (2004), the latter based on high 
signal to noise Space Telescope Imaging Spectroph (STIS) 
observations and extrapolated to higher wavelengths using 
Kurucz model atmospheres. After normalising the synthetic 
spectra to the SED of Vega we measured the B — V colours 
by applying the selected filter transmission curves. 

The results are shown in Fig. 1191 We also plot, for 
comparison, the photometric version of this colour following 
Vazdekis et al. (1996) which is computed with the Alonso et 
al. (1996, 1999) calibrations. The figure shows a good agree- 
ment for nearly all ages and metallicities, within the typi- 
cal photometric uncertainties derived from the zero point of 
the B and V filters, as shown by the colours differences us- 
ing two different normalizations of the Vega spectrum. We 
obtain rather similar time evolution colour curves for all 
the metallicities, with an offset of ~0.02 mag. For the SSPs 
with [M/H] — —2.32 the offset can be as large as ~0.05mag 
(MILES SEDs colours slightly bluer), which is well below 
the photometric errors of the colours. Note, however, that 
for this metallicity the MILES models are barely acceptable, 
as shown by the Q parameter (see Fig. [6|. This is because 
of the stellar parameter coverage of the MILES library and 
the fact that colour transformation are not accurate at such 
low metallicity. 

A rather similar offset (~0. 015 mag.) has been quoted 
in Paper I (see their Fig. 8) by comparing the B — V colours 
derived from MILES stellar spectra to the values tabulated 
in the Lausanne photometric database (Mermilliod et al. 
1997) for the same stars. As discussed in Paper I this off- 
set sets an upper limit to the systematic uncertainties of 
our photometry. For the sake of clarity we do not show in 
Fig. [I9]the colours obtained from the model SEDs computed 
with the temperature scale of Gonzalez-Hernandez & Boni- 
facio (2009). These models provide redder B — V colour by 
~0.02 mag, i.e. similar to the zero point uncertainties. This 
is because for a given temperature the stellar population 
model selects cooler spectra as the Gonzalez-Hernandez & 
Bonifacio (2009) temperatures are hotter than MILES by 
51 K. 

The results presented here show that the colours derived 
from our models, either those predicted via photometric li- 
braries or the ones measured from our MILES SSP SEDs, 
are consistent. Such agreement shows the reliability of the 
temperatures adopted in Paper II, which are consistent with 
the temperature scales of Alonso et al. (1996,1999), i.e. the 
photometric libraries that feed our models. We also provide, 
in electronic form and on the web (see Section [7]) , colours 
measured directly in the MILES spectra using combinations 
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Figure 19. We plot tiic broad band B — V colour derived from 
the SSP SEDs for different ages and metallicities (as indicated 
within the panel) and the Kroupa Universal IMF with the zero 
point set with two different Vega spectra in black (very thin) 
and red (thicker) (see the text for details). The thickest green 
line represents the photometric colour computed by Vazdekis et 
al. (1996), as updated in this work, which is based on extensive 
empirical photometric stellar libraries. 



of various filter systems that are commonly used in the lit- 
erature (e.g., SDSS, HST). 



6 APPLICATIONS TO STELLAR CLUSTER 
AND GALAXY DATA 

We compare, in this section, our new model SEDs with a set 
of representative stellar clusters of varying age and metallic- 
ity. We also illustrate the use and capabilities of these models 
for the study of galaxies with varying age, metallicity, abun- 
dance ratio and mass, by giving a few examples for some 
well-studied galaxies. 

6.1 Stellar Clusters 

Galactic globular clusters allow us to uniquely test our evolu- 
tionary synthesis models, since independent age and metal- 
licity estimates are available through deep CMD analyses. 
These studies have concluded that, in most cases, globular 
clusters can be considered as formed in a single burst, which 
is characterized by a single age and a single metallicity, i.e. 
an SSP. Note, however, that there are clusters for which two 
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Figure 20. Spectra of various representative Galactic globular 
clusters with varying metallicity, from Schiavon et al. (2005), plot- 
ted in black. For each cluster we overplot in green the best fitting 
SSP SED. We adopt for the models a Kroupa Universal IMF. The 
age and metallicity of the SSP are indicated within each panel. 
The SSP SEDs were smoothed to match the instrumental reso- 
lution of the observations (i.e. cr = 120kms~^). The obtained 
residuals are plotted at the same scale. Note that no continuum 
removal was applied when performing the fit. 



stellar components have been identified (e.g., Meylan 2003) 
and that stochastic effects (e.g., Cervino et al. 2002) may 
affect significantly their integrated spectra. 

We have selected three representative Galactic globu- 
lar clusters, with varying metallicity, from the sample of 
Schiavon et al. (2005), which is composed of 40 clusters. 
These authors obtained integrated spectra of high signal- 
to-noise (S/N--250) covering the range AA3350-6430 A with 
3.1 A (FHWM) resolution. Fig.EOlshows the spectra and the 
best SSP SED fits overplotted. We used, for this purpose, 
our SSP SED library with Kroupa Universal IMF. The SSP 
SEDs were conveniently smoothed to match the instrumen- 
tal resolution of the observations (i.e. a ~ 120 km s^^). We 
normalized the stellar cluster spectra and the model SEDs 
at 5000 A, but we did not modify the spectral shape of the 
data (and the models) to perform the fits. The best fit is 
understood as the one that minimizes the residuals, which 
are also plotted in the figure at the same scale. The obtained 
ages and metallicities fully agree with the CMD-derived es- 



timates of the Harris (1996) catalog. For 47Tuc we obtained 
the same fit as in V99. For NGC 5904 we derived a younger 
age, in agreement with the very deep CMD study of Marin- 
Franch et al. (2009) who have estimated for this cluster a 
relative age of 0.83, i.e. the cluster is located within their 
young GC branch. We also found a younger age for this clus- 
ter, by means of a line-strength analysis that included the 
newly defined H/3o age indicator of Cervantes & Vazdekis 
(2009). The higher H/? value obtained for this cluster is at- 
tributed in part to the presence of an increased fraction of 
Blue Stragglers (see Cenarro et al. 2008). 

The residuals of Fig. [20] do not show any significant 
colour difference in the covered spectral range, but some 
very low frequency variations that can be fully attributed 
to differences in the flux calibration. In fact, note that these 
low frequency residuals are qualitatively identical among the 
three clusters. However the residuals show differences in var- 
ious molecular bands, mainly in the blue part of the spectra 
for all the clusters. In particular, CN features at --^3860 A and 
~4150 A, are much stronger in the clusters than in the mod- 
els. We also see a residual in the G band at ~4300 A, being 
this absorption weaker in the data. The CN-strong absorp- 
tion, and its anticorrelation with CH, is a well known feature 
in the Galactic globular cluster family, which has been shown 
in both their individual stars (e.g., Norris & Freeman 1979) 
and its integrated spectra (e.g. Rose & Tripicco 1986). A sim- 
ilar result is obtained in V99 for 47Tuc, but using a different 
input stellar spectral library and also a different spectrum for 
the cluster. Such CN-strong absorption feature has also been 
reported for extragalactic globular clusters, like NGC 1407 
(Cenarro et al. 2007b). Their Lick/IDS C24668 vs. CN2 di- 
agram shows the same anticorrelation with respect to the 
scaled-solar abundance ratios. Overall, we find that the resid- 
uals blueward of ^4300A are significantly larger than for the 
redder part of the spectra for all the clusters. We do not 
discard possible effects due to the fact that the Milky Way 
globular clusters show oxygen-enhanced abundance ratios for 
low metallicities. In fact Cassisi et al. (2004) have shown that 
such effects become relevant for the spectral ranges covered 
by the B, and mostly U, broadband filters. Therefore an ap- 
propriate modelling for these clusters requires working with 
such a-enhanced stellar evolutionary isochrones as well as 
the use of stellar spectra with similar abundance pattern. In 
addition we also require to use specific stellar spectra with 
CN-strong absorption features and CN-enhanced isochrones. 

In Fig. [5T]we show our SED fit to the integrated spec- 
trum of the standard open cluster M 67 (Schiavon, Caldwell 
& Rose 2004). Unhke for the globular clusters of Fig. [20] 
these authors obtained the integrated spectrum for M 67 by 
co-adding individual spectra of cluster members, weighted 
according to their luminosities and relative numbers. The 
obtained age and metallicity is in good agreement with these 
authors (see also Schiavon 2007) and with our isochrone fit- 
ting results (Vazdekis et al. 1996) within the uncertainties. 
The residuals only show a low frequency pattern, which is 
related to differences in the fiux calibration. However the 
higher frequency residuals are negligible. It is particularly 
interesting that, unlike in Fig. 1201 we do not see the residu- 
als blueward of 4300 A. This shows that all relevant spectral 
types are present in the MILES library. Therefore the residu- 
als obtained for Galactic globular clusters can be attributed 
to the non-solar abundance ratios present in these clusters. 
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Figure 21. The spectrum of the Milky Way open cluster M67 
(Schiavon ot al. 2004) plotted in black. We ovcrplot the best fitting 
SSP SED in green, with its age and metallicity indicated within 
the panel. We adopt a Kroupa Universal IMF. The model has been 
smoothed to match the instrumental resolution of the observed 
spectrum (also quoted). The residuals are plotted in the bottom 
panel with the same scale. Note that no continuum removal was 
applied for performing the fit. 
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Figure 22. The spectrum a bright globular cluster of the galaxy 
NGC 1407, from Ccnarro ct al. (2007b), plotted in black. We over- 
plot the best fitting SSP SED in green, with its age and metallicity 
indicated within the panel. We adopt a Kroupa Universal IMF. 
The model has been smoothed to match the instrumental reso- 
lution of the observed spectrum (also quoted). The residuals are 
plotted in the bottom panel with the same scale. Note that no 
continuum removal was applied for performing the fit. 

However we propose to use our base models in the same way 
that we describe in Section r6.2.1l to deal with such abundance 
patterns. 

We illustrate in Fig. [22] the potential use of these 
models for the study of extragalactic globular cluster sys- 
tems. We plot a Keck/LRIS spectrum of a globular cluster 
of the massive elliptical galaxy NGC 1407 (Cenarro et al. 
2007b), with mean S/N~40. The model SEDs were properly 
smoothed to match the resolution of the observed spectrum 
{a lOOkms"^). We obtain an age of 17.78 Gyr and a 
metallicity of [M/H]= —0.71 in agreement, within the er- 
rors, with the values derived from the detailed line-strength 
analysis performed by these authors. 

6.1.1 Using line-strength indices to estimate GC ages and 
metallicities 

The line-strengths derived from our base models have al- 
ready been compared to the line-strengths of the integrated 



spectra of Galactic globular clusters (e.g., Mendel et al. 2007; 
Cenarro et al. 2008; Cervantes & Vazdekis 2009). These com- 
parisons confirm the well known model zero-point problem 
affecting the age estimates obtained from the Balmer lines, 
i.e. the derived ages are larger than the CMD-derived ages. 
Note that the same problem appears when using other mod- 
els. For an extensive description and discussion of this prob- 
lem we refer the reader to the above papers and to, e.g., 
Gibson et al. (1999), Vazdekis et al. (2001a), and Schiavon 
et al (2002). 

In Fig. [23] we plot several line-strength diagnostic dia- 
grams to obtain the ages and metallicities of the Galactic 
globular cluster sample of Schiavon et al. (2005). We use, for 
this purpose, our reference models and the alternative set of 
models described in Section [2. 3. II which adopt the recently 
published temperature scale of Gonzalez-Hernandez & Boni- 
facio (2009) . We see that the data fall outside the model grid 
towards ages larger than the oldest predicted SSP (see also 
Fig. I12|l . This is older than the ages derived from CMDs. If 
we use the model grid that adopt the Gonzalez-Hernandez 
& Bonifacio (2009) temperature scale we obtain, however, 
lower H/3o values, i.e., we alleviate in part the model zero- 
point problem obtaining slighly younger ages for the globu- 
lar clusters than with the reference models, using the MILEs 
temperatures. This is because when synthesizing a SSP SED 
the code assigns a cooler spectrum (i.e. with lower Balmer 
index values) to a requested star of given parameters, as the 
Gonzalez-Hernandez & Bonifacio (2009) temperatures are 
~51 K hotter than in MILES. Note also that these models 
provide slightly larger values for the metallicity indices. How- 
ever, as pointed out in Section [2]3TT1 we do not find a strong 
physical reason to use the Gonzalez-Hernandez & Bonifacio 
(2009) temperature scale, particularly for metal-rich objects. 
Note that Fig. [23] also shows that a significantly larger tem- 
perature offset (--^150 K), i.e., much larger than the current 
uncertainties in the temperature scale, is required to be able 
to overcome the model zero-point. Therefore other aspects, 
such as those discussed in Vazdekis et al. (2001a) and Schi- 
avon et al. (2002) should be considered. 

It is remarkable the large age sensitivity of the H/3o 
index, which allows us to obtain nearly orthogonal model 
grids, particulalry for [M/H]> —0.71. We see that our age 
estimates for the most metal-rich globular clusters do not 
depend on the metallicity indicator in use. However the ob- 
tained metallicities might differ among the plotted diagnos- 
tic diagrams, indicating a departure from the scaled-solar 
abundance pattern. In fact, based on the H/3o vs. CN2 dia- 
gram we obtain for many globular clusters metallicities that 
are significantly larger than the ones inferred from the H/3o 
vs. [MgFe] plot. Therefore these clusters show CN2 values 
that are larger than expected according to their total metal- 
licities, as traced by the [MgFe] index (e.g. Thomas et al. 
2003). Note that characteristic CN-strong pattern emerged 
in the residuals obtained from the full-spectrum fits shown 
in Fig. [50] 

Finally, note that these index-index diagrams show that 
there are a group of clusters with apparently younger ages. 
Cenarro et al. (2008) found a clear correspondence between 
the H/3 values and the specific frequency of Blue Stragglers 
(BSS), such that the higher the BSS ratio the younger the 
apparent spectroscopic age. We refer the reader to that paper 
for further details. 
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Figure 23. The H/3o age indicator of Cervantes & Vazdekis (2009) is plotted versus various mctallicity indices, all measured at LIS-5.0A 
resolution, for the Galactic globular cluster sample of Schiavon et al. (2005). Overplotted are our reference models (solid lines), which 
adopt the MILES temperatures, and another set of models that adopt the recently published temperature scale of Gonzalez-Hernandez 
& Bonifacio (2009) (dotted lines). SSP metallicities are quoted on the top of the model grid, just above the line corresponding to the 
youngest stellar populations (i.e. 3.16 Gyr). We indicate the age of the SSPs (in Gyr) for the most metal-rich models (i.e.[M/H]= -1-0.22). 
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Figure 24. Various Virgo elliptical galaxies from the sample of Yamada et al. (2006), with S/N >100 at the H7 feature, are plotted 
in black. For each galaxy we overplot in green the SSP SED, with Kroupa Universal IMF, corresponding to the age and metallicity 
(indicated within each panel), which have bee derived by these authors via the H^ja- vs. [MgFe] diagnostic diagram. The model spectra 
were smoothed to match the total resolution of the spectra (also indicated within the panels). Note that no continuum correction was 
applied. The residuals are plotted within the panels with the same scale. 



6.2 Galaxies 



In this section we show the use of our models for studying 
galaxy spectra. With this aim we have chosen four galaxies 
(out of 14) of the Virgo early-type galaxy sample of Yamada 
et al. (2006), with spectra of extraordinarily high S/N (S/N 
>100A~^ at the H7 feature), extracted within an aperture 
of rJlO. Our subsample of galaxies (NGC4387, NGC4473, 
NGC4489 and NGC4621) cover a range in age, metallicity, 
abundance ratio and mass, as these galaxies were selected 
alongside the Colour-Magnitue Relation of the Virgo clus- 



ter. The resolution of the Subaru spectrum for NGC4489 
is FWHM=3.lA, whereas for the other three galaxies, ob- 
served at the William Herschel Telescope (Vazdekis et al. 
2001b), the resolution is FWHM=2.4A. Such high quality 
spectra were obtained to be able to measure the Ji^a age in- 
dicator of Vazdekis & Arimoto (1999), which is virtually in- 
sensitive to metallicity. Mean luminosity weighted ages and 
metallicities were obtained by plotting £[70- versus the to- 
tal metallicity indicator [MgFe] (Gonzalez 1993; Vazdekis et 
al. 2001b; Thomas et al. 2003). The derived values are tabu- 
lated in Table 5 of Yamada et al. (2006). We did not attempt 



here to fit the spectra of these galaxies, but in Fig. 1241 we 
overplot for each galaxy the SSP SED with the closest age 
and metallicity, as listed in their Table 5. The model SEDs 
were smoothed to match the total resolution of the galaxy 
spectra (i.e. taking into account the instrumental resolution 
and galaxy velocity dispersion), which are indicated in the 
panels. 

The residuals of the comparison for the low- velocity dis- 
persion (a) , old galaxy, NGC 4387, are very small. For the 
other galaxy with low a but young age, NGC 4489, the 
residuals in the blue part are as small as those obtained 
for NGC 4387. However, we see a pattern in the residu- 
als at longer wavelengths. If the stars of this galaxy would 
have formed in two bursts then a characteristic residual pat- 
tern would be seen, namely a surplus of light in the blue 
and a deficit in the red, which is not the case. The good 
flux-calibration of the models presented here allows this 
type of analysis. Unlike the other galaxies, the spectrum of 
NGC 4489 was taken in one of the Subaru observing runs, 
for which the correction of the continuum shape to a fiux 
scale is less accurate than is achieved for the observing runs 
at the WHT. The same applies to the feature seen in the 
residuals at ~4225 A. 
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Figure 25. [ZMg/Zpc] vs. [Mg/Fc] abundance ratio estimates for 
the Virgo elliptical galaxies of Yamada et al. (2006). The proxy 
for the abundance ratio is the metallicity difference derived from 
plotting the H7o- age indicator of Vazdekis & Arimoto (1999) 
vs. the magnesium-dominated index, Mgb, and H'ja vs. the iron- 
dominated index, Fe3. We plot the [ZMg/^Fe] values listed in 
Table 5 of Yamada et al. To estimate the [Mg/Fe] abundance ra- 
tio we use the Lick/IDS indices and the models of Thomas et al. 
(2003). 



6.2.1 Dealing with abundance ratios using the base models 

For the most massive galaxies, NGC 4473 and NGC 4621 we 
see in Fig. [24] positive residuals for the Fe features, e.g., 
Fe4383, Fe5270, Fe5335, and negative residuals for the Mg 
feature at ~5170 A. This shows the well known [Mg/Fe] en- 
hancement, which is characteristic of massive galaxies (e.g., 
Peletier 1989; Worthey, Faber & Gonzalez 1992; Gonzalez 
1993). 

We show here how our base models can be safely used for 
studying galaxies with non-solar element mixtures. A good 
example of such use is shown in the line-strength analysis 
of Yamada et al. (2006). A proxy for the [Mg/Fe] abun- 
dance ratio is obtained with the aid of the base model grids 
and an age indicator that is almost insensitive to metallic- 
ity. If a galaxy spectrum is enhanced in Mg/Fe, we obtain a 
higher metallicity when this age indicator is plotted vs. a Mg- 
dominated index than when the index is Fe-dominated. Ya- 
mada et al. use the Jiy^ index of Vazdekis & Arimoto (1999), 
which requires spectra with extremely high S/N, and the 
metaUicity indicators Mg6 and Fe3 {= f°4383+Fc527o+Fc5335 
defined in Kuntschner, 2000). Irrespective of the metallic- 
ity indicator in use, such optimized age indicators provide 
orthogonal model grids, and the lines of constant age are es- 
sentially horizontal, which means that a given measurement 
for the age indicator corresponds to a unique age (see Cer- 
vantes & Vazdekis 2009). The metallicity difference obtained 
in this way, [ZMg/Zpe], is a good proxy for the abundance 
ratio [Mg/Fe], as was shown in, e.g., Sanchez-Blazquez et 
al. (2006b), de la Rosa et al. (2007) and Michielsen et al. 
(2008). In Fig.[5S]we plot the values obtained by Yamada et 
al. (2006, Table 5) for this proxy vs. the [Mg/Fe] estimates 
derived with the aid of the models of Thomas et al. (2003), 
which specifically account for non scaled-solar element par- 
titions. To apply these models we use the Lick/IDS version 
of the H^, Mg6 and Fe3 indices of Yamada et al. (2006). As 
an example. Table 5 of Yamada et al. lists [ZMg/Zpe]~0.6 for 
NGC 4473 and NGC 4621. Following Fig. [25] this translates 



to [Mg/Fe] ~ +0.25 on the abundance ratio scale of Thomas 
et al. (2003). Note that this [Mg/Fe] abundance scale is also 
a considerable approximation because i) the use of scaled- 
solar isochrones, ii) it is based on three synthetic spectra 
to model the SSPs and Hi) a fixed element ratio pattern is 
adopted. Note also that in Fig. [25] we combine the results 
from two different model predictions, and that our method 
requires extrapolating to higher metallicities to obtain Ziiig 
(see the corresponding plots in Yamada et al. 2006). This 
explains in part the larger dispersion obtained for the higher 
abundance ratios. 

It is worth noting that the newly defined H/3o index 
of Cervantes & Vazdekis (2009) provides similar age disen- 
tangling power as the 117^ index, but with the advantage 
that H;9o requires spectra of similar qualities as usually em- 
ployed for measuring the standard H/3 index (see Cervantes 
& Vazdekis 2009). H/3 can also be used for this purpose, but 
the resulting model grids are significantly less orthogonal, 
increasing the error bars of the [Zivig/Zpc] determinations. 

Finally, to obtain the total metallicity, these optimized 
age indicators can be plotted versus the combined indices 
[MgFe] (Gonzalez 1993) or [MgFe]' (Thomas et al. 2003), 
which have been shown to be virtually insensitive to the 
[Mg/Fe] abundance ratio. 

In the second paper of this series we will present an 
updated version of Fig. 1251 where we use both the scaled- 
solar and the a-enhanced version of our MILES models, in 
order to be fully consistent. In addition, these models reach 
higher metallicities, avoiding in this way the extrapolation 
that might be required to obtain [ZMg/Zpc]. However, we 
have shown here that our base models, which combine scaled- 
solar isochrones with an empirical stellar library, can be used 
for studying galaxies with a-enhanced element partitions via 
the line-strength analyses. 



6.2.2 Cautionary remark on the Hp model predictions 

Figure [26] shows the very popular diagnostic diagram H/? vs. 
[MgFe], which is obtained from our MILES SSP SEDs once 
smoothed to match the LIS-14.0A resolution. For compari- 
son we show the equivalent model grid from the V99 models, 
as updated in this work, which use the Jones (1999) stellar 
spectral library. Fig. [26] shows that the V99 model grid is 
significantly more orthogonal than the grid obtained with 
the MILES SSP SEDs, mainly due to the lower sensitivity 
of R/3 index to metallicity variations. The MILES SSP grid 
looks, in fact, more similar to the equivalent standard grid 
obtained on the basis of the empirical fitting functions of 
Worthey et al. (1994) on the Lick/IDS system. 

We plot the radial line-strength profile for the Virgo 
galaxy NGC 4387, which has been obtained from low resolu- 
tion Keck/LRIS spectra by Sanchez-Blazquez et al. (2007). 
We deliberately do not plot the observational errors, as we 
are only interested in illustrating the net effect on the trend. 
Note that we also plotted in Fig. [24] the higher resolution 
spectrum obtained at the William Herschel Telescope for 
the central aperture re/10, i.e., 1.6arcsec (Yamada et al. 
2006). According to the V99 model grid, the observed gra- 
dient can be attributed to a decreasing metallicity with in- 
creasing galactocentric distance. Using the MILES SSP grid, 
the gradient can be explained by a combination of increasing 
age and decreasing metallicity. 

We have investigated the possible cause for the differ- 
ences between the two model grids. As we use here an up- 
dated version for the V99 models, which is now based on 
the same isochrones and transformations to the observa- 
tional plane that we employ for the MILES models, these 
two ingredients cannot be the the ones we are looking for. 
We identified three other possible sources: i) intrinsic differ- 
ences in the spectra of MILES and Jones (1999) stellar li- 
braries, mostly around the H/3 feature (such differences have 
been quoted before by Worthey & Ottaviani (1997); ii) dif- 
ferences in the stellar parameters adopted here (i.e. Paper II) 
and in V99 and in ) differences in the algorithms employed for 
computing/assigning a stellar spectrum to a given set of at- 
mospheric parameters. In order to understand this issue, we 
synthetized SEDs using the 245 stars, both in the Jones and 
MILES libraries, with the following characteristics: (a) Mod- 
els using both stellar libraries. MILES and Jones. (b) Models 
using MILES spectra assigning them the two set of stellar 
parameters, those used in V99 and in MILES (c) Models us- 
ing MILES stellar spectra and with both algorthims, the one 
used in V99 and the new algortihm described in V03. 

For the sake of brevity we do not show here this detailed 
study but we summarize the main results. We found that 
these model differences are mostly originated by the stellar 
spectrum assignation scheme employed in V99, which works 
at a significantly more coarse stellar parameter resolution. 
The other aspects described above also have a significant 
impact, but to a lesser degree. Furthermore the II/3 index 
measurements on the SSP SEDs synthesized with the same 
algorithm but exchanging the stellar spectral libraries are 
consistent with the H/3 index measurements obtained on a 
star by star basis. Unfortunately, this problem might have 
some impact on certain studies that make use of the H/3 
index derived from V99 models. 

Given the improvement on the stellar spectrum assig- 
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Figure 26. H/3 vs. [MgFe] diagnostic diagram obtained from 
our MILES base model SEDs (solid line) and V99 (dotted line), 
which uses the Jones (1999) stellar library. The indices were mea- 
sured on the SSP spectra once smoothed to match the LIS-14.0A 
resolution. We plot the radial line-strength profile for the Virgo 
galaxy NGC 4387, which has been obtained from low resolution 
Keck/LRIS spectra by Sanchez-Blazquez et al. (2007). The size of 
the filled circles decrease with increasing galactocentric distance: 
1.25, 2.5, 5, 10, 20arcsec. 



nation scheme used for the MILES SSPs, we will provide an 
updated version of the V99 predictions in our website using 
the new algorithms. 



7 THE WEB TOOL 

The exploitation of stellar population synthesis models has 
depended in part on the development of tools by users. Some 
authors offer software to either help to interpret or aid the 
user to handle the provided predictions. A few examples 
include LECTOR (V99), GALAXEV (Bruzual & Chariot 
2003), PEGASE (Fioc & Rocca-Volmerange 1997,1999), or 
EZ Ages (Graves & Schiavon 2008). In an effort to continue 
and extend these initiatives, we have created a new website 
('http://miles.iac.es ) to provide support for using the spec- 
tra of both the MILES stellar library and the model predic- 
tions described in this paper, but also in V99 at 1.8 A resolu- 
tion. The webpage also handles the near-IR stellar spectral 
library CaT and models that we published in Cenarro et 
al. (2001a), Cenarro et al. (2001b), Cenarro et al. (2002) 
and V03. Line-strength predictions based on empirical fit- 
ting functions (including the ones of the Lick/IDS system) 
are also included. Furthermore the webpage includes photo- 
metric quantities such as fluxes, colours, mass-to-light ratios 
and surface brightness fluctuations, which have been pub- 
lished in Vazdekis et al. (1996) and Blakeslee, Vazdekis & 
Ajhar (2001). It is very important to note that all these 
predictions will be continuously updated with the latest de- 



velopments and with the feedback from users. These comply 
with the standards of the Virtual Observatory. 

This webpage is intended for users not only to retrieve 
both the stellar spectra and the synthetic SSP SED predic- 
tions in this project, but also to provide webtools to facilitate 
the handling and transformation of the spectra. Beyond the 
basic tools for the exploration and visualization of the whole 
dataset, at a first stage, we also provide advanced webtools 
to: 

• Generate synthetic stellar spectra for a given set of stel- 
lar atmosphere parameters. These spectra are computed on 
the basis of our stellar libraries according to the mathemati- 
cal algorithm described in Appendix B of V03 (see also Sec- 
tion [13]). 

• Transform the stellar spectra or the SSP spectra to the 
instrumental setup of the data observed by the user. This 
tool will allow the user to input the details of the instrumen- 
tal configuration of the observed data (spectral resolution 
and sampling) and generate a set of spectra matching that 
setup. The model/library spectra will be convolved with the 
observational PSF (typically Gaussian) to match the spec- 
tral resolution of the observed data. During this step, partic- 
ular attention will be paid as to whether the user requires the 
FWHM to be constant as a function of wavelength (in which 
case convolution will be performed in linear space), or con- 
stant in velocity (in which case convolution will be performed 
in logarithmic space). Finally, the data will be rebinned (i.e. 
spectrally sampled, and if desired, in log scale). 

• Measure a variety of line-strength indices and colors on 
the spectra. We make use of LECTOR (V99) and INDEXF 
(Cardiel 2007) for this purpose. This tool allows to input 
newly defined indices by the user. 

• Convert Lick/IDS line-strength indices to the new Line 
Index System proposed here. We use for this purpose the 
transformations described in Appendix A. 

• Estimate the best luminosity weighted age and metal- 
licity, and alpha-enhancement for a given set of line indices 
using RMODEL (Cardiel et al. 2003). This program deter- 
mines stellar population parameters (e.g. age, metallicity, 
IMF type and slope), using as input line-strength indices by 
interpolating SSP model predictions. Both linear and bivari- 
ate fits are computed to perform the interpolation. 

• Provide spectra for a given star formation history. 

These tools are intended to be the first of many more 
sophisticated tools, that we will incorporate in the website 
over time, to facilitate the comparison between the observed 
data and our stellar library and models, and in this way 
make the analysis easier. 



8 SUMMARY AND CONCLUSIONS 

Here we present stellar population model SEDs based on the 
empirical stellar spectral library MILES . In Paper I we have 
presented and described this library, whereas in Paper II we 
have published the stellar parameter determinations. Here 
we present synthetic SEDs for single-age, single-metallicity 
stellar populations (SSPs) covering the full optical spec- 
tral range, 3540.5-7409.6 A, at moderately high resolution, 
FWHM=2.3A, which is virtually constant as a function of 
wavelength (see Fig. 4 of Paper I). The SSP SEDs have a 



reliable, accurate flux calibration, which is one of the ma- 
jor advantages of the MILES library. Furthermore the red- 
der part of the SSP spectra do not show telluric residuals, 
as these were carefully cleaned from the stellar spectra as 
shown in Paper I. 

We consider the SEDs presented here as our base mod- 
els, as they combine an empirical library, which is im- 
printed with the chemical composition of the solar neigh- 
bourhood, with scaled-solar stellar isochrones (Girardi et 
al. 2000). Therefore, whereas our models are self-consistent, 
and scaled-solar for solar metallicity, this is not the case for 
the low metallicity regime, as the observed stars there do 
not show this abundance pattern (see Schiavon 2007). In 
a second paper we will present self-consistent models, both 
scaled-solar and a-enhanced, for a range in metallicities. For 
this purpose we have used MILES, together with theoretical 
model atmospheres, which are coupled to the appropriate 
stellar isochrones. 

An interesting feature of our base models is that they 
rely, as much as possible, on empirical ingredients. Apart 
from using MILES stellar spectra, we also use transforma- 
tions of the theoretical parameters of the isochrones to obser- 
vational, measurable, quantities that are based on extensive 
empirical, photometric libraries, with marginal dependence 
on model atmospheres. The implementation of MILES in 
the population synthesis code has been done very carefully: 
each star of the library has been compared with other MILES 
stars with similar atmospheric parameters, using the same 
interpolation algorithm used in the synthesis code. As a re- 
sult of these tests and other considerations, we have removed 
or decreased the weight with which a given star contributes 
to the synthesis of a stellar spectrum for 135 (out of 985) 
stars of MILES. 

The unprecedented stellar parameter coverage of the 
MILES library has allowed us to safely extend our optical 
SSP SEDs from intermediate to very old age regimes (0.06- 
18 Gyr) and the metallicity coverage of the SSPs from su- 
persolar ([M/H] = -1-0.22) to [M/H] = -2.32. The very low 
metallicity SSPs should be of particular interest for globular 
cluster studies. In addition, we have computed SSP SEDs 
for a suite of IMF shapes and slopes: unimodal and bimodal 
IMF shapes as defined in Vazdekis et al. (1996) with slopes 
varying in the range 0.3-3.3 (Salpeter: unimodal with slope 
1.3), and the segmented IMFs of Kroupa (2001) (Univer- 
sal and Revised cases). In addition to our reference models, 
which adopt the MILES temperatures published in Paper II, 
we provide an alternative set of models that adopt the re- 
cently published temperature scale of Gonzalez-Hernandez 
& Bonifacio (2009) (i.e., about ~51K hotter). The most rel- 
evant feature of these models is that the well known model 
zero-point problem affecting the age estimates obtained from 
the Balmer lines, i.e. the obtained ages are larger than the 
CMD-derived ages, is partially alleviated. 

To constrain the range of ages, metallicities and IMFs 
where the use of these SEDs are safe, we provide a quan- 
titative analysis. This analysis, which basically takes into 
account the parameter coverage of the stellar library feed- 
ing the models, clearly shows that our SSP SEDs can be 
safely used in the age range 0.06-18 Gyr, and for all the 
metallicities, except for [M/H] = —2.32, for which only stel- 
lar populations with ages above ~ 10 Gyr can be considered 
safe. This applies to the standard IMF shapes. However our 



quantitative analysis has shown that the quahty of the SSP 
SEDs decreases as we increase the IMF slope, especially for 
the Unimodal IMF. For example, for slopes above 1.8, the 
synthesized SSP SEDs for [M/H] = -2.32 are no longer safe 
for a Unimodal IMF. According to this analysis our pre- 
dictions are of significantly higher quality than any previ- 
ous model predictions in the literature for all the ages and 
metallicities that have been computed in this work. For ex- 
ample, any predictions based on the standard fitting func- 
tions of the Lick/IDS system (Worthey et al. 1994), e.g. line- 
strength indices, with which most stellar population studies 
have been performed so far, should not be considered safe for 
[M/H] < -1.3 (for standard IMF shapes). This limit is set 
at higher metallicity, [M/H] ~ —0.7, for predictions based 
on STELIB (e.g., Bruzual & Chariot 2003). The comparison 
of our new models, based on MILES, with the SSP SEDs 
published in V99, reveals systematic differences due to limi- 
tations in the flux calibration of the Jones (1999) stellar li- 
brary that feeds the latter models. As for the line-strengths 
the most important difference between these two models is 
found for the H/3 index, which shows a lower metallicity sen- 
sitivity in the V99 models. We show an example where this 
effect has a non-negligible impact on our interpretation of 
galaxy age /metallicity gradients. 

The SSP SEDs increase our ability to perform the stel- 
lar populations analysis, which can be done in a very flexi- 
ble manner. Observed spectra can be studied by fitting the 
full spectrum or by focussing on selected line indices, ei- 
ther standard or newly defined. An example of the previous 
method can be found in Koleva et al. (2008), where various 
stellar populations models were compared following this ap- 
proach. Our favourite approach is to analyze the spectrum 
of a galaxy at the resolution imposed by its velocity disper- 
sion (and its instrumental resolution), which require us to 
smooth the SSP SEDs to match the observed total resolu- 
tion. Therefore, the line index analysis to obtain relevant 
stellar population parameters is performed in the system of 
the galaxy, either through full spectrum fitting or line index 
approaches. 

As many authors prefer to plot the line-strength mea- 
surements of their whole galaxy sample in a single index- 
index diagnostic diagram, or to compare their data to tab- 
ulated indices in the literature, we propose a new Line In- 
dex System (LIS), which makes this method straightforward 
to apply. With LIS we avoid the well known uncertainties 
inherent to the widely employed Lick/IDS system. Unlike 
the latter, LIS has two main advantages, i.e., a constant 
resolution as function of wavelength and a universal flux- 
calibrated spectral response. Data can be analyzed in this 
system at three different resolutions: 5 A, 8.4 A and 14 A 
(FWHM), i.e. a =127, 214 and 357kms"^ at 5000 A, re- 
spectively. These resolutions are appropriate for studying 
globular cluster, low and intermediate-mass galaxies, and 
massive galaxies, respectively. Line-strength measurements 
given at a higher LIS resolution can be smoothed to match 
a lower LIS resolution. Furthermore we provide polynomials 
to transform current Lick/IDS line-index measurements in 
the literature to the new system. We provide LIS line-index 
tables for various popular samples of Galactic globular clus- 
ters and galaxies. We also show various popular index-index 
diagnostic diagrams for these samples in the LIS system. 

As an application, we have fitted a number of represen- 



tative stellar clusters of varying ages and metallicities with 
our models, obtaining good agreement with CMD determi- 
nations. Unlike for the open cluster M 67, our SED fits for 
a sample of Galactic globular clusters show non negligible 
residuals blueward of 4300 A, which mostly reflect the char- 
acteristic CN-strong features of these clusters among other 
deviations from the scaled-solar pattern. We also applied our 
models to representative galaxies with high quality spec- 
tra, for which independent studies are available, obtaining 
good results. We show that our base models can be used for 
studying line-strength indices of galaxies with a-enhanced 
element partitions. Examples of such a use can be found 
in, e.g., Vazdekis et al. (2001b), Kuntschner et al. (2002), 
Carretero et al. (2004), Yamada et al. (2006). The method 
consists in plotting a highly sensitive age indicator, such as 
those of Vazdekis & Arimoto (1999), or the recently defined 
H/3o index of Cervantes & Vazdekis (2009) versus Mg (e.g., 
Mg&) and versus Fe (e.g., Fe4383, (Fe)), obtaining a virtu- 
ally orthogonal model grid where the estimated age does not 
depend on the metallicity indicator in use. Unlike the age, 
for a [Mg/Fe] enhanced galaxy the obtained metallicities dif- 
fer when plotted against, for example, Mgfo and (Fe) indices. 
This metallicity difference, [ZMgb/Z<Fc>], can be used as a 
good proxy for the abundance ratio determined with the 
aid of stellar population models that specifically take into 
account non-solar element partitions. Although this proxy 
yields larger [Mg/Fe] values, there is a linear relation be- 
tween these two ways of estimating the abundance ratios 
(e.g., Sanchez-Blazquez et al. 2006b; de la Rosa et al. 2007; 
Michielsen et al. 2008). In the second paper of this series 
we will present both a-enhanced and scaled-solar self con- 
sistent SSP SEDs, all based on MILES, but modified with 
the aid of model atmospheres. By applying these models to 
galaxies with varying abundance ratios, for which high S/N 
spectra are available, we confirm the linear relation men- 
tioned above. There is one advantage in using this proxy: 
the abundance ratios determined in this way do not depend 
on the specific details of modeling of a-enhanced stellar pop- 
ulations such as those that are implicit to the computations 
of the atmosphere, or the adopted element mixtures. 

Our fits to the stellar cluster and galaxy spectra of high 
quality have confirmed the high precision of the flux calibra- 
tion of our SSP SEDs. Furthermore the colours that we ob- 
tained from our SSP SEDs are consistent with the ones that 
we derived via the employed photometric stellar libraries, 
within typical zero point uncertaintes. Such agreement shows 
the reliability of the temperatures adopted for the MILES 
stars (Paper II), which are consistent with the temperature 
scales of Alonso et al. (1996,1999). Such accuracy opens new 
applications for these models. For example, the model SEDs 
can be used as templates for codes aiming at determining 
spectrophotometric redshifts, particularly for those based on 
narrow-band filters (e.g. ALHAMBRA, PAU). 

Finally we present a webpage from which not only these 
models and stellar libraries can be downloaded, but where 
we also provide a suite of user-friendly webtools to facil- 
itate the handling and transformation of the spectra. For 
example, once the user enters the details of the instrumental 
setup employed in the observations (e.g., PSF, sampling), 
the users can obtain their favourite line strength indices 
and diagnostic diagrams, ready to plot their observational 
measurements. Apart of the model SEDs shown here, the 



webpage also provides predictions for a suite of observable 
quantities. 
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0.002 


C4668 


0.642 


0.831 


0.017 




H/3 


0.116 


1.001 


-0.003 


0.001 


Fe5015 


-0.154 


1.162 


-0.025 




Mgi 


-0.001 


0.886 


0.146 


0.140 


Mg2 


-0.004 


0.924 


0.013 


-0.028 


Mgb 


0.116 


0.950 


0.009 


-0.001 


Fe5270 


0.256 


0.705 


0.122 


-0.013 


Fc5335 


0.259 


0.805 


0.141 


-0.019 


Fg5406 


0.289 


0.629 


0.266 


-0.047 


Fe5709 


0.006 


1.045 


-0.052 


0.016 


Fe5782 


0.119 


0.467 


0.693 


-0.213 


Na5895 


-0.074 


0.921 


0.023 


-0.001 


TiOi 


-0.004 


0.773 


0.723 


-0.826 


Ti02 


0.001 


0.924 


0.116 


-0.132 


H<5a 


-0.369 


1.012 


1.012 




H7A 


-0.151 


0.974 


0.001 




HSf 


0.028 


0.970 


0.037 


-0.003 




-0.174 


1.050 


0.013 


-0.002 



Table Al. Polynomials to transform from the Liek/IDS system 
to the new system defined here for GC (LIS-5.0A) (see Eg. IA1|| . 



Ilis — ao + ailLick/iDS + a'2lLick/iDS + a-alLick/iDS (Al) 

where Ilis and iLick/iDS represent the index values in the 
LIS (at either of the three resolutions) and Lick/IDS system, 
respectively. Figure lXTI shows the results of the comparison 
between the indices of the Lick/IDS library and of MILES 
at FWHM=5A resolution. Similarly, Fig. |J2] and Fig. \M\ 
show the comparison of the indices of the Lick/IDS stars 
and of the MILES stars broadened to FWHM=8.4A and 
FWHM=14 A, respectively. The coefficients of these polino- 
mia, which convert the Lick/IDS indices to the LIS-5.0A, 
LIS-8.4A and LIS-14.oA, are given in Table Table \M\ 
and Table IA3I respectively. 



APPENDIX A: TRANSFORMATIONS FROM 
THE LICK/IDS SYSTEM TO THE NEW LINE 
INDEX SYSTEM 

In this Appendix we provide the required transformations 
that make it possible to convert Lick/IDS index measure- 
ments to the new system proposed here. For this purpose we 
have compared the index measurements of the 218 MILES 
stars in common with the Lick library. The indices were 
measured once the MILES stellar spectra were broadened 
to 5.0 A, 8.4 A and 14.0 A to obtain the conversions corre- 
sponding to LIS-5.0A, LIS-8.4A and LIS-14.0A, respectively. 
To calculate these conversions we fitted S"^"^ order polynomi- 
als in all cases. These polynomials are defined as follows: 




Hi^(Lick) 



H7p(Lick) 



Figure Al. Comparison of the Lick/IDS indices of the Lick/IDS library and those in MILES broadened to a constant resolution 
FWHM=5 A. Solid lines show the transformation from one system to the other. 




Figure A2. Comparison between the Lick/IDS indices of the Lick library and those measured on the stars in common with the MILES 
library broadened to FWHM=8.4A. Solid lines represent the transformation from one system to the other (see Tabic [AjJ 
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Figure A3. Comparison between the Lick/IDS indices of the Lick hbrary and those measured on the stars in common with the MILES 
library broadened to FWHM=14A. Solid lines represent the transformation from one system to the other (see Table [A3|l 
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TiOl 


-0.004 


0.766 


0, 


.736 


-0, 


.841 


TiOa 


0.001 


0.919 


0, 


.128 


-0, 


.142 


H<5a 


-0.246 


0.977 


0, 


.006 






H7A 


-0.175 


0.986 










H5f 


0.027 


0.973 


0, 


.016 


-0, 


.001 


H7f 


-0.149 


1.010 


0, 


.015 


-0, 


.002 



Table A2. Polynomials to transform from the Lick/IDS system 
to the new system defined here for dwarf and intermediate-mass 
galaxies (LIS-8.4A). 





ao 


ai 




a2 


as 


CNi 


0.005 


0.875 


-0, 


.091 


0.402 


CN2 


-0.002 


0.863 


-0, 


.110 


0.477 


Ca4227 


0.134 


0.380 


0, 


.113 


-0.008 


G4300 


0.251 


0.910 


0, 


.019 


-0.004 


Fe4383 


-0.203 


0.749 


-0, 


.004 


0.001 


Ca4455 


0.031 


0.331 


0, 


.073 


0.005 


Fe4531 


0.211 


0.814 


-0, 


.025 


0.004 


C4668 


0.500 


0.779 


0, 


.006 




H/3 


0.257 


0.861 


-0, 


.003 


0.001 


Fe5015 


0.123 


0.738 


-0, 


.004 




Mgi 


-0.003 


0.885 


-0, 


.164 


0.856 


Mg2 


-0.003 


0.874 


0, 


.154 


-0.134 


Mgb 


0.112 


0.738 


0, 


.022 


-0.001 


Fe5270 


0.224 


0.565 


0, 


.076 


-0.007 


Fe5335 


0.178 


0.496 


0, 


.090 


-0.010 


Fe5406 


0.160 


0.547 


0, 


.084 


-0.015 


Fe5709 


-0.018 


0.929 


-0, 


.064 


-0.013 


Fe5782 


0.038 


0.490 


0, 


.254 


-0.056 


Na5895 


-0.022 


0.801 


0, 


.034 


-0.002 


TiOi 


-0.004 


0.749 


0, 


.798 


-0.902 


TiOz 


0.000 


0.909 


0, 


.144 


-0.155 


Ha^ 


-0.017 


0.914 


0, 


.008 




H7A 


-0.126 


0.971 








B.Sf 


-0.074 


0.909 


0, 


.013 


-0.001 


H7f 


-0.199 


0.925 


0, 


.014 


-0.002 



Table A3. Polynomials to transform from the Lick/IDS system 
to the new system defined here for massive galaxies (LIS-14.0A). 



